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1 .  INTRODUCTION 


1.1  Summary 

Primary  activities  during  1981  consisted  of  (i)  expanding  the 
laboratory  capabilities  for  processing  and  testing  of  polymeric  compos¬ 
ites,  especially  by  integrating  automated  data  acquisition  equipment  with 
laboratory  instruments  and  developing  software  for  acquiring  and  analyz¬ 
ing  data;  (ii)  conducting  research  in  accordance  with  the  Statement  of 
Work  given  in  Section  1.2;  (iii)  preparing  nine  technical  papers  and 
M.S.  theses;  and  (iv)  various  interactions  of  the  faculty  with  the  tech¬ 
nical  community  through  presentation  of  papers,  participation  as  members 
of  technical  committees,  etc.  Reference  [1]  summarizes  work  done  in  1980. 

Sections  2-7  summarize  the  research  activities.  The  professional 
personnel  associated  with  the  project  and  outside  activities  of  the 
faculty  related  to  composites  are  given  in  Section  8.  Additional  material 
on  completed  work  and  student  activities  is  in  the  Appendices. 

1.2  Statement  of  Work 

"a.  Investigate  the  effects  of  cure  cycle  parameters  on  mechanical 
characteristics  of  resins,  composites  and  composite  structural 
specimens: 

(1)  Study  the  curing  process. 

(2)  Investigate  the  effects  of  cure  cycle  parameters  on  physical 
aging. 

(3)  Investigate  generation,  relief  and  effects  of  residual 
stresses  in  laminates. 

b.  Investigate  deformation,  damage  growth  and  fracture  behavior  for 
resins,  composites  and  composite  structural  specimens: 

(1)  Develop  and  verify  constitutive  equations. 

(2)  Develop  and  verify  ply  damage  and  delamination  models. 

(3)  Develop  an  automated  structural  material  property  characteri¬ 
zation  system  and  data  base." 
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2.  MATRIX  CONTROLLED  FRACTURE  BEHAVIOR 
OF  GRAPHITE/ EPOXY  COMPOSITES* 

2.1  Introduction 

It  is  now  recognized  that  delamination  cracking  is  a  serious  problem 
for  some  applications  of  graphite/epoxy  composites.  It  is  therefore  impor¬ 
tant  to  characterize  the  delamination  fracture  toughness.  Typically,  both 
tensile  and  shear  stresses  exist  on  the  fracture  plane,  and  thus  various 
combinations  of  tension  and  shear  have  to  be  studied  in  mixed-mode  frac¬ 
ture  tests.  Furthermore,  a  better  understanding  of  how  the  fracture  of 
the  epoxy  is  effected  by  the  presence  of  the  graphite  fibers  is  needed  so 
intrinsic  toughness  in  the  material  system  may  be  fully  realized  by 
proper  adjustments  in  processing  that  effect  the  thickness  of  the  resin 
rich  region,  fiber  waviness,  etc. 

The  purpose  of  this  research  has  been  to  better  understand  the  frac¬ 
ture  behavior  of  graphite/epoxy  composites  when  the  crack  path  is  princi¬ 
pally  through  the  epoxy  matrix.  Macroscopic  measurements  of  fracture 
toughness  have  been  made  to  determine  the  critical  energy  release  rate  for 
crack  propagation.  These  have  been  interpreted  in  terms  of  microscopic 
mechanisms  of  fracture  determined  from  a  careful  examination  of  the  frac¬ 
ture  surface  using  a  scanning  electron  microscope. 

The  majority  of  the  research  in  this  area  is  detailed  in  the  two  M.S. 
Theses  (cf.  Section  7  for  Abstracts): 

(i)  "Mode  I-Mode  II  Delamination  Fracture  Toughness  of  a  Unidirectional 
Graphite/ Epoxy  Composite,"  by  P.  S.  Vanderkley. 


♦Prepared  by  W.  L.  Bradley 
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(ii)  "Mode  I  Transverse  Cracking  in  Epoxy  and  a  Graphite/ Epoxy  Com¬ 
posite,"  by  D.  R.  Williams. 

2.2  Experimental  Procedure 

Two  types  of  specimens  have  been  employed  in  this  work;  namely,  com¬ 
pact  tension  specimens  (CT)  and  double  cantilever  beam  specimens  (DCB), 
as  shown  in  Figures  1  and  2,  respectively.  Panels  of  Hercules  3502, 

12"  x  6"  x  0.1",  have  been  made  with  and  without  AS1  carbon  fibers.  The 
slot  with  the  crack  starting  notch  was  made  parallel  to  the  fibers  for  the 
CT  specimens  cut  from  the  panel  containing  unidirectional  fibers.  Thus, 
the  crack  path  was  transverse  to  the  loading  direction  but  parallel  to  the 
fiber  direction;  i.e.,  as  in  transverse  cracking  in  a  90°  specimen.  The 
CT  specimens  of  neat  material  and  graphite/epoxy  were  precracked  by  gently 
tapping  a  razor  blade  placed  at  the  root  of  the  machined-in  notch,  inducing 
a  natural  crack  to  form. 

Similar  graphite/epoxy  panels  12"  x  6"  x  0.1"  were  prepared  with  a 
strip  of  teflon  placed  along  one  edge  at  the  midthickness.  Subsequently, 

1"  x  12"  x  0.1"  DCB  specimens  were  cut  from  these  panels  with  the  teflon 
tape  serving  as  a  crack  starter  at  one  end.  When  these  specimens  are 
loaded  symmetrically,  pure  Mode  I  delamination  results.  However,  when  they 
are  asymmetrically  loaded,  a  combined  mixed  mode  I/II  fracture  results,  as 
shown  in  Figure  3. 

It  should  be  noted  that  the  transverse  cracking  in  the  CT  specimen  is 
not  identical  to  the  delamination  cracking  in  the  DCB  specimen,  though  both 
pass  principally  through  the  matrix.  The  delamination  cracking  occurs 
through  a  resin  rich  region  which  occurs  between  layers  and  is  typically 
1  to  2  fiber  diameters  in  thickness.  The  transverse  cracking  parallel  to 
the  fibers  occurs  through  whatever  resin  is  available  in  the  respective 


Figure  1.  Scale  drawing  of  compact  tension  specimens  used 
in  this  study,  with  dimensions  in  mm. 


Load  Cell 


specimen,  and  grip  configuration 


Figure  3.  Superposition  of  a  combined  loaded  beam  into  a  symmetric 
(pure  mode  I),  and  an  asymmetric  (pure  mode  1 1)" case 
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plies,  though  this  is  considerably  less  than  that  available  locally  between 
plies  for  the  usual  method  of  processing  laminates.  This  difference  will 
be  seen  to  significantly  affect  the  fracture  toughness  results  for  matrix 
cracking. 


2.3  Analysis  of  Experimental  Results 

The  experimental  results  for  the  asymmetrically  loaded  DCB  specimens 
have  been  analyzed  using  linear  elastic  fracture  mechanics  coupled  with 
linear  beam  theory,  which  gives  [2]: 


GI  = 


2  2 
8PsLc 
BEI 


’ll 


6PaLc 

~wr 


(2.1) 


where  B  is  the  DCB  width,  El  is  the  stiffness  of  the  beam  section  that  is 

not  yet  cracked,  L_  is  the  crack  length  in  the  DCB,  P  is  the  symmetric 

P  +  P  P  -  P. 

Toad,  — and  Pa  is  the  asymmetric  Toad,  —  ,  where  and  PL  are 

the  upper  and  lower  loads  as  seen  in  Figure  3.  The  total  energy  release 

rate  is  (Gj  +  Gjj).  Note  that  when  the  loading  is  symmetric  (i.e., 

Pu  =  PL)  Pa  .  0  and  G„  .  0. 

The  critical  energy  release  rate  was  determined  in  two  ways  from  the 
compact  tension  specimens.  First,  compliance  measurements  were  made  for 
specimens  with  various  relative  crack  lengths,  a/w,  and  then  the  critical 
load  to  fracture  the  specimen  was  determined.  The  critical  energy  release 
rate  was  found  using  the  relationship  [3] 

Glc'^§  <2-2> 


where  Pg  is  the  critical  load,  B  is  the  thickness  of  the  compact  tension 
specimen  and  C  is  the  specimen  compliance  (5/P),  where  5  =  the  total  load 
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line  displacement  and  P  is  the  applied  load. 

The  second  method  of  determining  the  critical  energy  release  rate 
from  the  CT  specimens  was  indirectly  through  a  determination  of  the  criti¬ 
cal  stress  intensity  factor,  K^.  This  was  done  using  relationships  for 
K-jc  as  a  function  of  Pg  and  (a/w)  by  Srawley  [4]  for  isotropic  materials 
and  Mandell  et  al.  [5]  for  othotropic  materials.  The  two  approaches  give 
very  consistent  results. 

Finally,  the  critical  stress  intensities  were  converted  into  critical 
energy  release  rates  in  the  usual  way  with  [3] 


G,  =  Klc(1  ~ 

I  r 


V2) 


(2.3) 


where  E*  for  the  isotropic  neat  material  is  the  Young's  modulus  E,  and 
for  the  orthotropic  material  is  calculated  from  a  compliance  tensor  using 
the  relationship  for  effective  modulus  [6]: 


E* 


bl  1^22  1/2  jj^22y 72  +  2bl 2^ b66  j  ^ 


(2.4) 


Material  constants  from  Tsai  and  Hahn  [7]  for  AS1/3501  (which  is  very 

2 

similar  to  AS1/3502)  were  used  to  evaluate  E*/(l  -v  ). 


2.4  Experimental  Results  and  Discussion 

The  critical  energy  release  rate  results  for  pure  mode  I  cracking  may 
be  summarized  as  follows: 

2 

A.  neat  material  G^c  =  69  N-m/m 

B.  delamination  cracking  2 

parallel  to  fibers  G^c  =  150  N-m/m 

2 

C.  transverse  cracking  G^c  =  225  N-m/m 
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These  results  clearly  indicate  that  the  presence  of  fibers  in  a  brittle 
matrix  enhances  the  toughness  in  a  significant  way.  Earlier  results  by 
Bascom  et  al . [8]  show  that  fiber  reinforcement  in  a  ductile  matrix  decreased 
the  delamination  toughness.  These  different  effects  can  be  rationalized 
in  the  following  way.  The  deformation  or  damage  zone  ahead  of  the  crack 
tip  in  a  ductile  or  tough  material  involves  a  significant  volume  of  highly 
deformed  material,  contributing  to  the  energy  dissipation  that  accompanies 
fracture.  Fiber  reinforcement  seems  to  reduce  the  size  of  this  "plastic" 
zone,  giving  a  net  decrease  in  the  fracture  toughness. 

In  a  weak,  brittle  material,  reinforcement  by  strong  brittle  fibers 
enhances  the  toughness.  The  "plastic"  zone  in  epoxy  3502  is  negligibly 
small,  so  that  there  is  no  reduction  in  "plastic"  zone  size  associated  with 
the  introduction  of  fibers.  However,  where  cracking  occurs  through  a  thin 
layer  o'  resin  bounded  by  fibers  with  even  a  modest  amount  of  waviness, 
some  interaction  of  the  crack  tip  with  the  fibers  is  unavoidable.  If  the 
fibers  are  much  stronger  than  the  matrix,  as  they  are  for  graphite/epoxy 
systems,  then  this  interaction  makes  crack  extension  more  difficult. 

Where  this  occurs,  either  fiber  breakage  or  multiple  initiation  results, 
the  latter  giving  rise  to  a  scalloped  appearance  of  the  epoxy  on  the  frac¬ 
tured  surface.  The  greater  G-jc  for  transverse  cracking  as  compared  to 
delamination  (225  N/m  vs  150  N/m)  could  be  interpreted  as  resulting  from 
the  thinner  layer  of  resin  through  which  the  cracking  occurs,  which  results 
in  more  frequent  interaction  with  the  fibers. 

For  mixed  mode  delamination,  typical  results  are  seen  in  Fig.  4. 

These  results  indicate  that  increasing  mode  II  loading  increases  the  total 
energy  release  rate  required  for  fracture.  The  shear  stress  acting  on  the 
delamination  plane  causes  the  plane  of  the  maximum  normal  stress  to  be~ 
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rotateu  relative  to  the  plane  of  the  lamina,  the  magnitude  of  this  rotation 
increasing  with  increasing  mode  II  shear.  Crack  extension  in  the  brittle 
epoxy  apparently  will  seek  to  follow  the  maximum  normal  stress  plane. 
However,  cracking  in  this  direction  may  only  proceed  a  short  distance  be¬ 
fore  the  crack  path  is  interrupted  by  the  presence  of  fibers.  Again  either 
fiber  breakage  or  multiple  crack  initiation  results,  giving  a  scalloped 
fracture  surface  [9].  The  increasing  Gj  with  increasing  mode  II  shear 
is  the  result  of  the  crack  continually  being  redirected  into  the  fibers. 
While  macroscopically  we  may  say  the  fracture  is  mixed  mode,  microscopic¬ 
ally  the  fracture  behavior  is  probably  that  of  a  whole  series  of  short 
pure  mode  I  fractures  on  the  maximum  normal  stress  plane  with  tearing 
to  coalesce  these  cracks,  giving  a  crack  path  that  in  two  dimensions  looks 
like  a  saw  blade  for  a  handsaw.  It  is  easy  to  see  how  increasing  mode  II 
changes  the  p^iod  of  the  scalloping,  giving  a  greater  total  energy  dissi¬ 
pation  in  fracture.  The  above  micromechanistic  model  is  strongly  supported 
by  fractographic  observations  made  in  a  scanning  electron  microscope  as 
presented  in  the  theses  of  Vanderkley  [2]  and  Williams  [9], 

In  summary,  matrix  dominated  cracking  is  greatly  influenced  by  the 
presence  of  fibers,  the  greater  the  crack  tip  interaction  with  the  fibers, 
the  greater  the  toughness.  Thinner  resin  rich  regions  and/or  mode  II  shear 
serve  to  enhance  this  interaction,  requiring  larger  values  of  energy  re¬ 
lease  rate  for  crack  extension.  It  should  be  noted  the  above  generaliza¬ 
tions  are  strictly  applicable  to  a  brittle  matrix  reinforced  with  stronger, 
brittle  fibers.  Different  behavior  would  be  expected  for  a  composite  with 
a  ductile  matrix. 
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3.  MATRIX  CONTROLLED  FRACTURE  ANALYSIS 
OF  FIBROUS  COMPOSITES* 

3.1  Introduction 

Experimental  results  described  in  Section  2  and  in  earlier  reports 
on  this  project  [1,  10]  are  helping  to  provide  the  basis  for  develop¬ 
ment  of  theoretical  models  of  damage  growth  and  ^racture  of  polymeric 
composite  materials.  In  this  section  we  discuss  additional  experimental 
work  and  theoretical  models  which  are  concerned  with  the  same  objective. 
The  following  three  M.S.  Theses,  whose  Abstracts  are  collected  in  Sec¬ 
tion  7,  cover  a  large  part  of  the  work  discussed  in  this  section: 

(i)  "Compression  Induced  Delamination  in  a  Unidirectional 
Graphite/ Epoxy  Composite,"  by  J.  W.  Earley. 

(ii)  "Mode  I  Delamination  of  Unidirectional  Graphite/ Epoxy  Com¬ 
posite  under  Complex  Load  Histories,"  by  J.  S.  Cullen. 

(iii)  "Analysis  of  the  Effect  of  Matrix  Degradation  on  Fatigue 

Behavior  of  a  Graphite/ Epoxy  Laminate,"  by  R.  T.  Arenburg. 

The  next  three  subsections  briefly  discuss  the  theses  in  the  order 
listed.  The  last  subsection,  3.5,  abstracts  a  recent  report  of  work  on 
this  project  which  is  concerned  with  a  theory  of  crack  growth  in  non¬ 
linear  viscoelastic  media.  Both  nonlinear  behavior  and  viscoelasticity 
have  been  taken  into  account  so  that  the  theory  will  be  general  enough 
to  apply  to  damage  growth  and  fracture  of  composites  exposed  to  a  wide 
range  of  service  environments  over  extended  periods  of  time. 

3.2  Compression  Induced  Delamination 

The  objective  of  the  experimental  phase  of  this  research  was  to 
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use  in-plane  compressive  loading  to  induce  delamination  in  Mode  I  frac¬ 
ture.  A  beam-column  specimen  geometry  was  employed,  as  shown  in  Fig. 

5. 

All  of  the  beams  were  fabricated  from  Hercules  AS! -3502  graphite/ 
epoxy.  The  standard  15.24  cm  wide  prepreg  tape  was  cut  into  30.48  cm 
lengths,  and  these  plies  were  used  to  fabricate  a  twelve-ply  unidirec¬ 
tional  (0°)  laminate.  One  strip  of  teflon  film  was  placed  at  the  mid¬ 
plane  to  create  the  initial  delamination.  The  panels  were  press  cured 
using  a  press  which  was  automatically  controlled  by  a  microprocessor, 
programmed  with  temperature,  pressure  and  vacuum  profiles.  The  panels 
were  cut  into  2.54  cm  by  30.48  cm  strips  and  side  and  end  fittings  were 
attached  to  each  specimen.  After  the  fittings  were  attached,  the  col¬ 
umn  was  pre-cracked  to  assure  an  initial  symmetric  crack  and  to  allevi¬ 
ate  any  artificial  effects  arising  from  the  method  of  inducing  the 
delamination. 

Two  other  specimens  were  fabricated  to  determine  independently  the 
critical  energy  release  rate  from  laminated  beams  split  at  one  end;  cf. 
Fig.  3B  and  [IT],  The  beams  were  fabricated  by  cutting  the  beam  columns 
in  half. 

A  constant  side  loading  Ps,  indicated  in  Fig.  5,  was  applied  and 

then  the  beam  column  was  compressed  in  an  Instron  test  machine  using  a 

-4 

constant  crosshead  rate  of  2  x  10  cm/sec.  The  specimens  were  tested 
in  a  dry  state  at  room  temperature. 

The  side  loading  provided  a  well-defined  out-of-plane  deflection 
and  served  to  stabilize  the  beam  against  global  buckling.  We  wanted  to 
check  the  basic  theory  with  a  minimum  of  extraneous  specimen-related 
factors;  therefore,  it  was  decided  to  introduce  out-of-plane  deflections 


L 


Side  Fitting 
(for  measuring  lateral 
deflection ) 


Grips  for  applying 
Load  P. 


End  Fittings  ,  mount  to  lower 
compression  apparatus 


Test  Specimen 
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in  this  manner  instead  of  depending  on  (random)  effects  of  initial 
imperfections  such  as  slight  initial  beam  curvature  and  load  misalign¬ 
ment. 


Figures  6  and  7  show  typical  results  for  axial  deflection  and 
energy  release  rate;  the  latter  was  predicted  after  using  lateral 
displacement  measurements  to  calculate  crack  length.  The  theory  uti¬ 
lized  the  assumptions  of  classical  beam  theory  and  accounted  for  the 
geometric  nonlinearity  due  to  large  rotations. 

Primary  conclusions  reached  from  the  study  by  J.  W.  Earley  are: 

1)  For  AS1-3502  graphite/epoxy  in  a  dry  state  at  room  tempera¬ 
ture,  classical  beam-column  analysis  linked  with  linear 
elastic  fracture  mechanics  is  applicable  for  determining 
axial  and  lateral  deflections,  strain  energy,  and  energy 
release  rate.  (Similar  agreement  for  more  complex  geome¬ 
tries  has  been  reported  by  others  [12].) 

2)  The  average  critical  energy  release  rate  obtained  from 
beam-columns  was  137  N-m/nr  (.78  in-lb/ in2),  which  is 
close  to  the  value  of  127  N-m/m2  found  from  end-split 
laminates. 

3)  The  energy  release  rate  during  crack  progagation  was  fairly 
constant  throughout  each  test. 

4)  Increased  lateral  load  did  not  have  any  significant  effect 
on  energy  release  rate,  although  it  did  increase  the  average 
crack  speed. 

5)  The  overall  fracture  path  generally  propagated  through  the 
resin  between  fibers,  although  failure  along  the  interface 
with  some  fiber  breakage  was  not  uncommon. 

6)  Hackles,  overlapping  platelets  of  epoxy,  and  cleavage-like 
micro-cracks  which  branch  near  fiber  epoxy  boundaries  were 
observed.  These  features  could  be  used  to  determine  the 
direction  of  damage  progression  in  composite  specimens. 


3.3  Delamination  Under  Complex  Load  Histories 

The  specimen  design  chosen  for  this  study  was  a  twelve  ply  2.54  cm 
x  25.4  cm  (1  in  x  10  in)  Hercules  AS1/3502  graphite/epoxy  coupon  speci- 


calculated  axial  deflection 
actual  axial  deflection 


Axial  Deflection  vs.  Axial  Load  for  Specimen  43 
(crack  propagation  portion  only,  test  direction 
right  to  left). 
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men  with  fibers  parallel  to  the  beam  axis.  During  lay-up,  teflon  was 
inserted  at  the  mid-plane  of  the  laminate  at  one  end  of  the  specimen  to 
provide  an  initial  delamination  crack.  Aluminum  tabs  were  bonded  to 
the  same  end  in  order  to  load  the  samples  as  shown  in  Fig.  3B;  the  left 
end  of  the  beam  was  unsupported  in  these  tests.  For  the  elevated 
temperature/humidity  testing,  the  desired  moisture  content  was  obtained 
by  placing  specimens  in  an  environmental  chamber  at  a  humidity  of  95%  RH 
and  temperature  of  70°C.  The  constant  humidity  was  maintained  by  a 
saturated  solution  of  potassium  sulfate  in  the  bottom  of  the  closed 
chamber.  Equilibrium  moisture  content  was  obtained  by  conditioning 
the  samples  for  six  weeks  prior  to  testing. 

The  specimens  were  tested  at  ambient  (55%  RH,  23.9°C)  and  elevated 
(95%  RH,  93.3°C)  environments  under  the  transient  and  cyclic  crack  open¬ 
ing  displacements  (COD)  A  U$  in  Fig.  3B)  illustrated  in  Figs.  8  and 
9.  An  MTS  servohydraul ic  testing  machine  interfaced  with  a  Digital 
Equipment  MINC-11  data  acquisition  laboratory  computer  system  was  used 
for  the  experiments.  An  environmental  chamber  maintained  the  desired 
conditions  during  testing. 

Figure  10  shows  typical  cyclic  test  data;  crack  growth  occurs 
along  the  upper  load  curve,  which  decreases  with  increasing  displace¬ 
ment.  Observe  that  upon  reloading,  a  larger  load  is  needed  to  initiate 
crack  growth  than  existed  at  the  start  of  unloading.  Crack  growth  was 
not  continuous  along  the  upper  curve,  but  instead  propagated  in  a 
stop-start  fashion.  This  stop-start  behavior  also  existed  in  the  con¬ 
stant  displacement  rate  tests,  as  may  be  seen  in  Fig.  11.  This  type  of 
unstable  crack  growth  has  been  reported  by  others  in  studies  of  resin 
fracture  [13-15]. 


ure  10.  Load  vs.  Crosshead  Displacement  for  a  5  Cycle  Test 
(Specimen  28). 
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Figure  11.  Crack  Length  vs.  Time  for  a 

Ramp  Test  (Specimen  25,  93.3°C 
and  Specimen  6,  23.3°C). 
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A  quasi-static  viscoelastic  analysis  reported  last  year  [1]  pre¬ 
dicted  a  curve  of  G  versus  crack  speed  a  which  has  a  negative  slope 
except  at  very  low  speeds.  Figure  12  illustrates  such  a  curve,  which 
may  help  to  explain  the  observed  unstable  crack  growth.  When  G  exceeds 
the  value  required  for  growth  at  point  "a",  G-j q ,  the  crack  propagates 
dynamically.  However,  for  quite  slowly  applied  displacements,  the 
available  G  decreases  rapidly  until  it  reaches  the  value  G-|a,  the 
"arrest  energy".  Unless  the  applied  load  is  removed,  the  crack  does  not 
truly  arrest  in  this  model,  but  rather  has  the  speed  at  point  "d".  The 
proposed  G  vs  a  behavior  is  indicated  by  the  arrows  in  Figure  12,  and 
is  typical  of  the  behavior  found  from  all  of  the  tests;  Figure  13  shows 
such  behavior. 

Further  analysis  of  the  data  is  needed  to  determine  if  the  crack 
actually  stops  propagating  during  this  cycle  of  "stop-start"  propagation, 
or  if  it  just  slows  to  the  point  "d"  in  Fig.  12  before  repeating  the 
cycle.  A  more  detailed  study  will  have  important  implications  for  the 
understanding  of  crack  tip  deformation  processes  and  development  of 
crack  growth  theory.  A  theory  based  on  cycles  of  true  initiation  and 
arrest  will  be  quite  different  (in  detail,  at  least)  than  a  theory  based 
on  continuous,  but  uneven,  growth. 

General  conclusions  that  may  be  drawn  from  J.  S.  Cullen's  thesis 

are: 

1)  The  maximum  ("initiation")  energy  G-j g  increases  with  increasing 
temperature  and  moisture  and  decreasing  applied  displacement 
rates. 

2)  The  stop-start  type  of  unstable  crack  growth  is  most  apparent 
at  the  highest  temperature  and  moisture  levels. 

3)  The  minimum  ("arrest")  energy  is  approximately  90"  of  G-i q  for  the 
composite  materials  and  test  conditions  employed.  (Differences 
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Figure  13.  Energy  Release  Rate  vs.  Crack  Speed  Showing 
Local  Unstable  Crack  Growth  (dG/da  <  0). 
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between  initiation  and  arrest  values  for  resins  may  be 
much  larger  [13-15]. 


3.4  Analysis  of  Effect  of  Matrix  Degradation  on  Deformation  and  Strength 


Tensile  failure  of  laminates  is  often  a  direct  result  of  delamina¬ 
tion  which  initiates  at  free  edges,  such  as  illustrated  in  Figure  14. 

R.  T.  Arenburg's  thesis  is  concerned  with  the  prediction  of  the  mechanical 
state  and  energy  release  rate  for  delamination  in  tensile  coupons  with 
various  amounts  of  delamination  and  matrix  degradation  (e.g.,  due  to  micro¬ 
cracking).  A  detailed  study  of  a  [±45/902^  graphite/epoxy  laminate  was 
made  using  theory  developed  in  the  thesis  and  experimental  data  reported 
elsewhere  [16].  Guided  by  both  theory  and  experiment,  symmetrical  delami¬ 
nations  between  the  45°  and  90°  plies  were  assumed  to  be  the  primary 
macroscopic  cracks,  and  therefore  energy  release  rate  in  only  this  mode 
of  fracture  was  analyzed. 

Figure  15  shows  the  effect  of  reduced  values  of  the  resin's  Young's 
modulus  Em  on  the  predicted  axial  modulus  of  the  laminate;  EREp  is  the 
latter  modulus  based  on  the  unreduced  matrix  modulus  Ev^...^.  Experi¬ 
mental  values  of  the  axial  mean  strain  significantly  increased  during 
fatigue  loading  at  constant  load  amplitude  and  R-value  [16];  typical 
behavior  is  sketched  in  Fig.  16,  in  which  the  rapid  increase  in  strain 
near  failure  is  due  to  delamination.  Viscoelastic  effects  without  matrix 


damage  were  far  too  small  to  explain  the  strain  increase  (even  well  before 
the  high  strain  rate  period  near  failure)  [16],  and  there  was  no  signifi¬ 
cant  amount  of  fiber  failure  until  the  final  stages  of  fatigue  failure. 
Therefore,  prior  to  delamination  the  increase  in  strain  was  assumed  to  be 
due  to  matrix  damage,  idealized  for  purposes  of  theoretical  analysis  by  a 
cycle  and  stress-dependent  matrix  modulus.  Using  the  finite  element 
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Figure  16.  Idealized  mean  strain  curves  for  a  [±45/90?] 
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results  for  modulus  in  Fig.  15  and  the  strain  data  in  Fig.  16,  it  was 
estimated  that  Em/ i a  1  1135  to  ™  t^ie  ran9e  fr0m  0*1  to  0.25  just 
before  the  rapid  upturn  in  the  strain  curves  if  the  predicted  and  experi¬ 
mental  strains  are  to  agree  at  this  point;  the  largest  amount  of  matrix 
damage  near  failure  (ratio  of  0.1)  corresponds  to  the  lowest  stress  (50% 
UTS). 


Figure  17  gives  the  theoretical  prediction  of  energy  release  rate  as 
a  function  of  crack  length/ply  thickness  for  three  different  values  of 
matrix  modulus;  "O'Brien's  method"  is  based  on  classical  lamination  theory. 
In  order  to  test  the  proposed  model  of  fatigue  failure,  a  representative 
critical  energy  release  rate  of  Gc  =  175  N/m  was  used  (along  with  the  esti¬ 
mated  amount  of  matrix  degradation)  to  predict  the  initial  flaw  size  that 
would  result  in  unstable  delamination;  values  close  to  the  ply  thickness 
were  calculated,  which  is  a  very  encouraging  result.  Furthermore,  with¬ 
out  matrix  degradation,  the  energy  release  rate  in  Fig.  17  is  far  below 
that  needed  to  cause  unstable  delamination,  regardless  of  the  size  of 
the  initial  delamination,  for  the  stress  levels  used  in  [16]. 

Several  conclusions  can  be  made  from  R.  T.  Arenburg's  analysis: 

1)  Matrix  degradation  modeled  by  a  reduced  matrix  modulus  through 
micromechanics  principles  shows  considerable  promise  as  a  means 
to  analyze  fatigue  behavior.  However,  the  study  conducted  was 
of  a  preliminary  nature  in  that  only  mean  stress  and  mean  strain 
responses  were  considered,  and  the  influence  of  stress  gradients 
on  matrix  damage  was  neglected. 

2)  In  evaluating  laminate  stiffness,  it  was  found  that  as  the  matrix 
modulus  was  reduced,  considerable  shear  deformations  occurred 
which  classical  lamination  theory  was  unable  to  model.  The  shear 
deformations,  which  are  associated  with  free  edge  effects,  were 
found  to  penetrate  a  considerable  distance  in  from  the  free  edge. 
For  the  smallest  value  of  the  matrix  modulus,  the  finite  element 
method  showed  that  the  free  edge  effects  were  extended  five 
laminate  thicknesses  into  the  interior  of  a  [±45/902]  laminate. 
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Figure  17.  Total  normalized  strain  energy  release 
rate  as  a  function  of  delamination 
length/ ply  thickness  for  a  constant 
stress  condition. 
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3)  A  parametric  study  of  free  edge  delaminations  as  a  function  of 
matrix  damage  was  performed  for  constant  maximum  strain  and 
constant  maximum  stress  fatigue  tests.  For  the  constant  strain 
fatigue  test  the  available  energy  for  delamination  was  predicted 
to  be  greatest  at  the  beginning  of  the  test  and  then  was  subse¬ 
quently  reduced  by  the  matrix  damage.  In  a  constant  stress 
fatigue  test,  the  effect  of  matrix  damage  was  to  increase  the 
available  energy  for  delamination.  As  illustrated  in  Fig.  18 
for  a  laminate  with  constant  critical  energy  release  rate  Gc, 
with  sufficient  matrix  damage  (decreasing  matrix  modulus)  a 
crack  of  initial  length  aQ  will  become  unstable. 

3.5  Generalized  Models  for  Damage  Growth  and  Fracture 

The  preceding  discussion  illustrated  the  importance  of  accounting 
for  distributed  damage  in  the  prediction  of  deformation  and  fracture  of 
composites.  Considerable  extension  beyond  the  use  of  a  single  parameter 
reflecting  the  damage,  such  as  a  uniform  reduced  matrix  modulus,  is 
necessary  to  achieve  a  sufficiently  general  theory  for  engineering  uses. 
Microcracking  (and  other  damage  mechanisms)  in  general  will  be  nonuni - 
formly  distributed,  produce  directional  softening,  and  depend  on  histories 
of  the  mechanical  and  physical  environments. 

As  an  aid  to  understanding  and  prediction  of  complex  damage  phenomena, 
a  viscoelastic  crack  growth  theory  for  nonlinear  media  was  recently  devel¬ 
oped  on  this  project  [17].  Starting  with  a  nonlinear  viscoelastic  consti¬ 
tutive  equation  in  the  form  of  a  modified  superposition  integral,  corre- 
•  spondence  principles  are  established  in  [17]  which  enable  viscoelastic 
stresses  and  deformations  to  be  easily  predicted  from  a  nonlinear  elas¬ 
ticity  solution;  one  correpsondence  principle  applies  with  stationary  or 
growing  cracks,  while  another  may  be  used  with  time-dependent  joining 
and  healing  of  initially  separate  surfaces.  The  first  of  these  corre¬ 
spondence  principles  leads  to  a  generalized  J-integral  for  use  in  pre¬ 
dicting  fracture  initiation  time  and  crack  speed.  An  important  step  in 
this  analysis  is  the  derivation  of  nonlinear,  non-singular  solutions  for 
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Figure  18.  Delamination  growth  criterion  for  a  constant 
maximum  stress  fatigue  test. 
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the  mechanical  state  of  the  neighborhood  of  the  crack  tip;  singular  solu¬ 
tions  result  in  physically  unacceptable  viscoelastic  fracture  criteria. 
Several  explicit  results  for  the  initiation  and  growth  rate  of  cracks  are 
obtained  in  [17]  assuming  the  strain  energy  density  for  the  associated 
elastic  material  is  a  homogeneous  function  of  the  strains  (i.e.  the  stress 
is  a  power-law  function  of  strain). 

Results  from  this  theory  have  been  used  in  a  study  of  the  relation 
between  crack  speed  and  the  J-integral  [18],  and  in  the  development  of 
damage  models  for  composites  with  emphasis  on  effects  of  particle  rein¬ 
forcement  [19]. 

Currently,  we  are  working  on  extending  the  damage  theory  to  laminates. 


\ 
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4.  STRESS  AND  DEFORMATION  ANALYSIS  OF 
RESIN  AND  COMPOSITE  MATERIALS* 

4.1  Residual  Stresses 

An  investigation  of  the  effects  of  moisture  and  temperature  on  resid¬ 
ual  stresses  in  composite  laminates  is  continuing.  A  major  portion  of  the 
effort  during  the  past  year  was  conducted  as  part  of  a  Ph.D.  dissertation 
project  "On  the  Effects  of  Post  Cure  Cool  Down  and  Environmental  Condition¬ 
ing  on  Residual  Stresses  in  Composite  Laminates"  by  Mr.  B.  D.  Harper. 

Solution  forms  which  relate  curvatures  in  unsymmetric  cross-ply  lami¬ 
nates  to  temperature  excursions  and  moisture  content,  including  viscoelas¬ 
tic  effects,  were  formulated  and  programmed.  Progress  was  made  concerning 
the  specimen  design,  regarding  both  plates  and  beam-like  coupons.  The 
design  strives  to  keep  the  thickness  dimension  small  enough  to  ensure  both 
measurable  curvatures  and  reasonably  short  moisture  saturation  times.  On 
the  other  hand  it  is  desirable  for  deflections  to  be  small  enough  to  per¬ 
mit  usage  of  linear  beam  and  plate  theory.  Several  lay-ups  were  considered 
with  staggered  0°  and  90°  plies  to  avoict  premature  cracking  due  to  resid¬ 
ual  stresses.  The  range  of  validity  of  the  classical  linear  plate  theory 
was  examined  by  analyzing  some  different  large  deformation  problems. 

An  improved  device  to  measure  curvatures  was  constructed  and  several 
environmental  chambers  were  assembled.  Also,  as  part  of  the  dissertation 
project,  preliminary  work  was  done  on  thj  curing  of  bi -material  coupons 
made  of  Hercules  3501-6  resin  and  aluminum  backing  in  preparation  for 
investigating  cure-shrinkage  effects  in  AS/3501-6  graphite/epoxy  composites. 


*Prepared  by  Y.  Weitsman 
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Some  earlier  work  on  moisture  effects  is  covered  in  the  recent  publi¬ 
cation  by  Y.  Weitsman,  "A  Rapidly  Convergent  Scheme  to  Compute  Moisture 
Profiles  in  Composite  Materials  under  Fluctuating  Ambient  Conditions,"  J. 
Composite  Materials,  Vol .  15  (July  1981),  pp.  349-358.  A  reprint  is 
included  in  Appendix  A. 

4.2  Optimization  of  Cooling  Histories 

A  theoretical  investigation  was  conducted  on  optimal  cool-down  paths 
for  minimizing  residual  stresses  in  cross-ply  laminates  and  adhesive  joints. 

report  by  Y.  Weitsman  and  B.  D.  Harper,  "Optimal  Cooling  of  Cross-Ply 
Composite  Laminates  and  Adhesive  Joints,"  describes  this  work,  and  is 
included  in  Appendix  A. 

The  investigation  was  based  upon  a  previous  formulation  of  the  opti¬ 
mal  cool-down  idea,  but  extended  the  scheme  to  include  the  technically 
important  case  of  temperature-dependent  transverse  and  longitudinal  coef¬ 
ficients  of  thermal  expansion,  and  a^.  Specific  calculations  were  con¬ 
ducted  for  epoxy/aluminum  joints  and  two  kinds  of  graphite/epoxy  composites 
using  the  best  currently  available  data.  It  was  noted  that  while  the  opti¬ 
mal  time-temperature  paths  did  not  vary  appreciably  for  temperature-dependent 
ay  and  aL,  a  further  reduction  was  found  in  the  residual  stresses  at  the 
termination  of  cool-down  when  the  specific  variations  of  thermal  expansion 
coefficients  were  considered. 

4.3  Shear  Deformation  in  Laminated  Plates 

An  investigation  of  shear  deformation  effects  is  reported  by  Mr.  B. 

A.  Coulter  in  the  M.S.  thesis  "Shear  Deformation  Effects  in  Composite 
Laminates".  This  investigation,  which  is  abstracted  in  Section  7,  includes 
a  study  of  elastic  effects  in  plates  subjected  to  several  types  of  normal 
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loadings;  it  was  shown  that  sharp  spatial  variations  of  external  loads 
may  cause  very  significant  shear  effects,  sufficient  to  invalidate  the 
results  of  classical  plate  theory.  Such  departures  from  classical  theory 
may  occur  even  for  thin  laminates.  In  addition,  a  viscoelastic  version 
of  the  shear  deformation  theory  was  employed  to  predict  material  damping 
due  to  the  time-dependent  response  of  the  resin  component.  Again,  it  was 
found  that  under  sharp  spatial  variations  in  external  loadings,  as  well 
as  for  thick  laminates,  the  classical  plate  theory  may  grossly  under-predict 
the  amount  of  material  damping  available  due  to  the  viscoelasticity  of  the 
resin. 

4.4  Analysis  of  Crazing 

An  inquiry  was  initiated  to  represent  damage  growth  in  polymeric 
resins  due  to  crazing  prior  to  cracking.  It  appears  that  in  moist  environ¬ 
ments  crazes  tend  to  extend  and  stretch  under  constant  loads, ..a  behavior 
that  is  probably  attributable  to  the  accentuated  viscoelasticity  of  the 
crazed  material  which  diminishes  its  load  carrying  capacity.  On  the  other 
hand  dry  crazes  ("air  crazes")  appear  to  grow  by  drawing  additional  mate¬ 
rial  from  the  polymerized  regions  across  sharply  defined  interfaces.  In 
this  case  the  energetics  of  the  phenomenon  must  account  for  energy  required 
to  affect  such  phase  transformation.  The  mechanics  formulation  of  these 
behaviors  requires  the  solution  of  some  well-defined  boundary  value  prob¬ 
lems,  which  are  currently  being  set-up. 
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5.  VISCOELASTIC  FINITE  ELEMENT  PROGRAM* 


The  nonlinear  viscoelastic  theory  of  Schapery  [20]  was  examined  in 
detail  during  1979  and  1980  as  a  possible  means  of  accounting  for  heredi¬ 
tary  behavior  within  the  context  of  an  incremental,  nonlinear  finite  ele¬ 
ment  program.  The  integral  equations  were  recast  into  incremental  form 
and  a  time-marching,  recursive  numerical  procedure  was  developed  for  evalu¬ 
ating  certain  time-dependent  strain  functions.  During  1980  and  1981,  this 
numerical  procedure  was  incorporated  into  the  nonlinear  finite  element  pro¬ 
gram,  AGGIE,  in  a  general  three-dimensional  format.  Several  applications 
were  made  to  uniaxial  and  biaxial  creep  problems  and  comparisons  were  made 
to  experiment  and  other  creep  models.  Results  of  some  of  these  comparisons 
are  reported  [21,  22]. 

The  viscoelastic  constitutive  model  for  isothermal  conditions  and 
infinitesmal  strains  is  given  by  [20], 
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and  A1  ',  8  and  tr  are  constants.  The  functions  a^,  a^  and  occj  are 
material  properties  and  are  in  general  functions  of  the  three  stress  invari 
ants.  The  quantity  is  a  reduced  time  variable  while  are  engineering 
strain  components  and  accj  are  components  of  the  Cauchy  stress  tensor.  The 
function  aJU)  is  the  transient  component  of  the  so  called  linear  viscoe¬ 
lastic  shear  creep  compliance.  The  first  term  in  Eq.  (5.1)  represents 
the  instantaneous  elastic  response  of  the  material. 

To  cast  Eq.  (5.1)  into  a  form  suitable  for  use  in  an  incremental 
finite  element  analysis,  its  time  derivative  is  taken,  and  Eq.  (5.2a)  is 
used  to  obtain. 
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Equation  (5.3)  is  rewritten  as 
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A  time  marching  procedure  for  evaluating  Xcd  and  Ycd  is  given  in  [21]. 

The  evaluation  of  the  two  integrals  Xcd  and  Ycd  given  by  Eqs.  (5.5) 
and  (5.6)  represents  a  key  element  in  evaluating  the  creep  strain  rate. 
Once  the  values  of  these  integrals  are  known,  the  creep  strain  rate  is 
evaluated  from  the  sum  of  the  second  and  third  terms  of  Eq.  (5.4)  as: 
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The  material  parameter  ag,  ag.  3cd,  and  aD  are  evaluated  from  stan¬ 
dard  creep  and  recovery  tests.  Details  are  provided  in  [20,  22].  In 
some  special  cases  it  is  possible  to  "back-out"  these  material  properties 
from  creep  equations  as  outlined  in  [21]. 

In  the  computer  program,  these  material  properties  may  be  input 
either  in  functional  form  or  piece-wise  linear,  tabular  form. 

In  the  finite  element  program,  the  creep  strains  predicted  by  Eq. 
(5.8)  are  superimposed  with  elastic  and  time-independent  plastic  strains 
predicted  by  an  incremental  plasticity  theory  using  a  kinematic-isotropic 
hardening  model  and  von  Mises  yield  condition.  Limited  comparison  with 
experimental  data  indicates  relatively  good  correlation  for  non  metals 
and  some  metals  at  elevated  temperature.  Final  numerical  results  are 
pending  on  the  completion  of  Sanders  work  [22]. 
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6.  STRUCTURE- PROPERTY  STUDIES  OF  RESIN* 

6.1  Differential  Scanning  Calorimeter  (DSC)  Investigations 

Several  other  workers  have  studied  physical  aging  of  thermoplastics 
(in  contrast  to  thermosets,  which  we  have  investigated)  using  the  DSC. 

They  find  that  the  shape  and  position  of  glass  transition  changes  as  the 
sample  is  aged.  Qualitatively,  this  is  understood  by  the  free  volume 
gradually  annealing  from  the  sample.  In  some  cases  they  found  quantita¬ 
tive  effects  that  they  could  relate  to  mechanical  properties.  We  endeav¬ 
ored  to  extend  their  type  of  study  to  a  graded  series  of  epoxy  resins. 
However,  the  attempt  to  measure  aging  of  a  series  of  epoxy  resins  by  means 
of  changes  in  the  specific  heat  near  the  glass  transition  temperature,  T  , 
has  been  abandoned,  at  least  temporarily.  The  effects  we  have  seen  are 
small  and  buried  in  experimental  noise  so  that  interpretation  has  been 
difficult.  This  result  parallels  the  small  effect  of  aging  observed  in 
another  part  of  this  program  in  the  mechanical  characteristics  of  one  of 
the  structural  epoxy  resins,  as  discussed  in  the  Third  Annual  Report  [1], 
The  change  from  thermoplastics  to  thermosets  has  two  deleterious 
effects  on  the  study.  First,  the  thermosets  we  studied  have  higher  glass 
transition  temperatures  than  the  thermoplastics  studied  in  the  literature. 
There  is  a  decrease  in  the  magnitude  of  the  change  of  specific  heat  as  T^ 
increases.  Secondly,  the  sharpness  of  the  transition  is  reduced;  this  may 
be  a  result  of  either  an  inhomogeneous  sample  or  be  inherent  in  all  thermo 
sets.  In  the  particular  series  of  epoxy  resins  we  used,  the  curing  agent 
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is  di ami  nodi  phenyl  methane,  which  is  a  high  melting  solid.  In  order  to 
mix  this  with  the  resin,  the  combination  must  be  heated  to  a  moderately 
high  temperature.  While  the  mixture  appears  to  be  well  mixed  and  visibly 
homogeneous,  there  can  be  heterogeneities  on  a  molecular  scale  that  would 
not  be  visible  to  the  eye.  Theoretically,  one  would  expect  the  equilibrium 
condition  for  an  uncured  mixture  to  be  homogeneous;  but  phase  separation 
is  quite  possible  as  the  cure  proceeds  since  even  small  positive  enthal¬ 
pies  of  interaction  can  become  important  with  curing  underway.  If  this 
phase  separation  were  to  occur  late  in  the  crosslinking  step,  then  the 
fineness  of  the  phases  would  result  in  a  material  that  is  apparently,  but 
not  actually,  homogeneous. 

There  is  considerable  discussion  in  the  literature  concerning  the 
general  presence  of  heterogeneities  in  both  thermosets  and  thermoplastics. 
Some  electron  microscope  images  show  a  texture  that  has  been  interpreted 
as  a  fundamental  heterogeneity  in  a  wide  variety  of  material  including 
epoxies.  Since  optical  aberrations  in  the  electron  microscope  can  cause 
a  simitar  texture  and  since  there  are  not  simple  experimental  tests  that 
will  distinguish  whether  a  texture  comes  from  these  aberrations  or  is 
physical,  the  question  is  still  open. 

In  any  case,  it  appears  as  though  the  change  in  specific  heats  is 
small  and  occurs  over  a  wider  temperature  range  in  our  samples  when  com¬ 
pared  to  the  thermoplastics  studied  in  the  literature.  These  factors  make 
the  small  differences  due  to  aging  even  harder  to  distinguish  from  back¬ 
ground  noise. 

In  addition,  the  epoxy  resins  are  sensitive  to  traces  of  moisture 
and  this  modifies  the  DSC  trace  as  much  as  aging.  We  endeavored  to  con¬ 
trol  the  moisture  environment  of  our  samples  but  some  anamolous  runs  are 


attributed  to  a  failure  to  control  moisture  adequately. 


The  combination  of  a  small  and  diffuse  effect  in  the  specific  heats 
and  the  inherent  scatter  in  the  DSC  measurements  have  kept  us  from  observ¬ 
ing  any  systematic  effects  of  aging  in  our  samples. 

Attempts  to  improve  the  precision  of  the  DSC  have  not  been  very  fruit 
ful .  Although  there  are  some  reports  in  the  literature  as  to  how  this  has 
been  done,  we  have  not  succeeded  in  obtaining  results  which  permit  the 
study  we  set  out  to  do.  Further  studies  of  physical  aging  will  be  post¬ 
poned  until  we  acquire  testing  equipment  which  is  sufficiently  sensitive 
to  changes  in  mechanical  or  thermal  properties. 

6.2  Toughened  Resins 

The  sunnier  of  1981  was  spent  at  the  Air  Force  Materials  Laboratory 
(AFML),  Wright  Aeronautical  Laboratories,  where  a  study  on  toughening 
mechanisms  in  high  temperature  resins  was  continued.  A  report  on  the  1980 
summer's  work  has  been  written,  accepted,  and  proofed  for  a  technical 
report,  "Toughening  of  ATX  Resins";  it  is  now  in  press. 

The  AFML  has  a  program  to  develop  a  generation  of  high  temperature 
resins.  A  wide  variety  of  resins  is  being  developed  in  which  rather  short 
chain  compounds  are  crossl inked  through  terminal  acetylene  groups.  One 
can  utilize  a  wide  variety  of  backbone  chains  with  the  same  crosslinking 
reaction,  can  achieve  high  thermal  stability,  and  can  start  with  rela¬ 
tively  low  molecular  weights.  Therefore  the  understanding  of  the  general 
properties  of  these  compounds  should  enable  selection  of  a  good  resin  to 
be  used  in  composites. 

These  materials  can  be  made  in  many  different  compositions,  but  tend 
to  be  quite  expensive  to  prepare  in  experimental  quantities.  It  is  there¬ 
fore  of  importance  to  understand  as  much  about  these  compounds  as  possible 
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before  initiating  a  large  scale  program.  The  flexibility  in  composition 
should  enable  choosing  a  very  attractive  resin  if  the  various  tradeoffs 
are  clearly  understood.  One  of  the  factors  limiting  the  success  of  such 
a  resin  is  the  toughness  of  the  resin  and  of  the  composite  made  from  it. 

The  relationship  of  these  quantities  from  fracture  behavior  to  the  molec¬ 
ular  structure  is  mostly  conjecture  at  this  time.  Yet  the  payoff  of  any 
understanding  is  great.  Therefore,  the  summer  of  1981  was  spent  at  AFML 
primarily  seeking  molecular  parameters  likely  to  be  related  to  toughness 
in  fractionated  samples  of  the  corresponding  thermoplastic  materials. 

Fractions  were  obtained  from  P-1700  and  P-1800  polysulfone  thermo¬ 
plastics,  which  are  similar  in  backbone  structure  to  the  ATS  resins,  and 
fractions  were  obtained  from  PPQ-401 ,  a  polypheny! qui noxal ine  thenmoplastic 
which  is  similar  to  the  backbone  structure  of  another  family  of  acetylene 
terminated  compounds  considered  as  a  potential  high  temperature  resin. 

Dilute  solution  viscosity  was  used  to  establish  the  random  coil  con¬ 
figuration  of  the  molecules  in  solution. 

The  viscosity  of  more  concentrated  solutions  showed  an  interaction 
at  very  low  apparent  entanglement  lengths.  It  was  observed  that  this  crit¬ 
ical  concentration  for  interaction  between  chains  varied  as  M"^  rather 
than  as  M-^  found  for  flexible  chains.  A  mathematical  model  for  this 
interaction  has  been  developed  and  is  in  preparation  for  publication.  The 
feature  responsible  for  this  behavior  in  the  critical  concentration  is, 
according  to  this  model,  the  rigidity  of  the  chains.  It  is  not  necessary 
for  a  rigid  chain  to  have  the  conformation  of  a  rigid  rod,  but  it  can  have 
the  conformation  of  a  random  walk  if  the  rotational  transitions  of  the 
bonds  is  very  slow. 

High  temperature  measurements  on  the  undiluted  polymers  indicate 
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that  at  these  high  temperature,  the  chains  are  flexible  and  entangle  as 
flexible  chains. 

The  relationship  of  these  observations  to  toughening  of  the  resins 
and  composites  has  yet  to  be  made;  but  the  conjecture  is  that  the  cross¬ 
link  spacing  must  be  no  less  than  the  distance  between  entanglements  in 
order  to  obtain  the  maximum  elongation  at  break.  Higher  crosslink  density 
would  make  the  material  more  brittle  and  produce  a  higher  glass  transition 
temperature.  Lower  crosslink  density  would  give  little  benefit. 

Several  papers  on  these  topics  are  in  various  stages  of  preparation 
as  part  of  my  effort  on  the  AFOSR  project. 

6.3  Viscoelastic  Behavior 

A  molecular  model  for  creep  is  in  the  early  stages  of  development. 

Rather  general  and  plausible  assumptions  can  lead  to  a  creep  function  of 
1  /3 

Aexpkt  '  ,  which  matches  a  variety  of  experimental  studies.  Other  expo¬ 
nents  are  possible  depending  upon  the  assumptions.  It  is  not  clear  whether 
this  type  of  creep  behavior  is  due  to  statistically  homogeneous  molecular 
motion  or  if  it  is  a  property  of  cracking,  crack  healing,  or  other  phe¬ 
nomena  on  a  scale  large  compared  to  molecular  dimensions. 


7.  GRADUATE  RESEARCH  ASSISTANT  ACTIVITIES 


7.1  Summary 

The  third  group  of  graduate  engineering  students  to  participate  on 
the  M.S.  level  in  the  AFOSR  research  project  entered  the  program  in 
September  1980.  They  now  have  completed  all  of  their  requirements  for  a 
Master  of  Science  degree,  and  the  results  of  the  research  are  reported 
in  the  following  theses: 

1.  Arenburg,  R.T.,  "Analysis  of  the  Effect  of  Matrix  Degradation 
on  Fatigue  Behavior  of  a  Graphite/Epoxy  Laminate." 

2.  Braswell,  J.L.,  "The  Effect  of  Geometry  on  the  Design  of 
Filament-Wound  Fiberglass  Tension  Lugs." 

3.  Coulter,  B.A.,  "Transverse  Shear  Deformation  Effects  in  Composite 
Laminated  Plates." 

4.  Cullen,  J.S.,  "Mode  I  Delamination  of  a  Unidirectional  Graphite/ 
Epoxy  Composite  Under  Complex  Load  Histories." 

5.  Earley,  J.W.,  "Compression  Induced  Delamination  in  a  Unidirec¬ 
tional  Graphite/Epoxy  Composite." 

6.  VanderKley,  P.S.,  "Mode  I  -  Mode  II  Delamination  Fracture  Tough¬ 
ness  of  a  Unidirectional  Graphite/ Epoxy  Composite." 

7.  Williams,  D.R.,  "Mode  I  Transverse  Cracking  in  an  Epoxy  and  a 
Graphite  Fiber  Reinforced  Epoxy." 

The  abstracts  are  collected  in  subsection  7.2.  Copies  of  the  theses 
will  be  provided  upon  written  request  to  the  Principal  Investigator  (R.A. 
Schapery) . 

The  current  group  consists  of  three  M.S.  and  four  Ph.D.  students. 
Besides  conducting  research  on  this  project,  the  students  are  involved  in 
academic  courses,  as  noted  in  the  program  announcement  in  Appendix  B. 


7.2  Abstracts  of  M.S.  Theses 
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ABSTRACT 

Analysis  of  the  Effect  of  Matrix  Degradation  on  Fatigue 
Behavior  of  a  Graphite/ Epoxy  Laminate.  (May  1982) 

Robert  Thomas  Arenburg,  B.C.E.T.,  Southern  Technical  Institute 
Chairman  of  Advisory  Committee:  Dr.  R.  A.  Schapery 

An  analytical  study  of  the  effect  of  fatigue  on  a  [+45/902^ 
graphite/epoxy  laminate  is  conducted  considering  two  levels  of  damage. 
The  first  one  is  microscopic  damage,  which  is  viewed  as  a  matrix  soft¬ 
ening  mechanism  that  is  accounted  for  by  the  use  of  micromechanics. 

The  second  level  of  damage  is  delamination.  An  evaluation  of  lamina¬ 
tion  theory  versus  a  finite  element  elasticity  solution  which  accounts 
for  free  edge  stresses  is  made  first.  Lamination  theory  is  found  to 
be  inadequate  due  to  its  inability  to  account  for  large  amounts  of 
shear  deformation  as  the  matrix  degrades.  Experimental  data  are  used 
to  evaluate  the  matrix  degradation  as  a  function  of  fatigue  cycles, 
and  then  a  parametric  study  of  the  strain  energy  release  rate  as  a 
function  of  delamination  size  at  various  degrees  of  matrix  degradation 
is  conducted.  For  a  constant  maximum  stress  fatigue  test,  the 
matrix  degradation  leads  to  an  increase  in  strain  energy  which  is  pre¬ 
dicted  to  eventually  cause  rapid  delamination.  Reasonably  good  agree¬ 
ment  between  theory  and  experimental  data  on  the  [±45/9023s  laminate 
is  demonstrated. 
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ABSTRACT 

The  Effect  of  Geometry  on  the  Design  of  Filament-Wound 
Fiberglass  Tension  Lugs.  (December  1981) 

James  Lee  Braswell,  Jr.,  M.S.,  Texas  A&M  University 
Chairman  of  Advisory  Coircnittee:  Dr.  Richard  M.  Alexander 

Finite  element  analysis  is  employed  to  characterize  the  mechanical 
behavior  of  filament-wound  fiberglass  racetrack  style  tension  lugs  as  a 
function  of  lug  geometry.  The  effect  on  ultimate  load  carrying  capa¬ 
bility  and  the  location  of  failure  origination  due  to  static  loads 
applied  uniaxially  to  S2/E773  glass/epoxy  specimens  is  considered.  Three 
configurations  are  studied,  each  having  identical  cross-sectional  areas. 
They  differ  only  in  the  path  the  fibers  traverse  around  the  pin.  In  the 
baseline  case,  the  path  is  semi-circular.  The  other  two  follow 
semi -elliptical  paths  with  the  major  diameter  of  the  ellipse  aligned 
with  the  direction  of  the  applied  load.  The  major  to  minor  diameter 
ratios  are  1.5  and  2.0.  Stress  distributions  at  various  locations  of 
interest  are  calculated  and  presented  graphically.  The  lug  failure  load 
and  fracture  initiation  site  is  predicted  for  each  configuration.  The 
predictions  are  based  on  the  maximum  stress  criterion  of  failure  and 
published  material  property  data.  From  these  analytical  results  it  is 
concluded  that  increased  structural  efficiency  can  be  obtained  by  taking 
advantage  of  the  shape  latitude  available  through  filament  winding. 

Experimental  data  used  to  verify  the  analysis  is  obtained  by  loading 
seven  specimens  of  each  lug  shape  to  failure.  Two  of  the  seven  specimens 
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are  strain  gaged  in  locations  compatible  with  the  finite  element  grid  to 
aid  in  checking  the  analytical  results.  Wei  bull  location  parameters  are 
determined  for  each  lug  shape  and  expected  failure  loads  are  compared  to 
those  predicted  by  the  analysis.  Reasonably  good  predictions  are  made 
for  two  of  the  three  lug  shapes  regarding  failure  strength  and  location 
of  failure  initiation.  Accurate  results  for  the  third  require  additional 
considerations  not  treated  here. 
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ABSTRACT 

Tranverse  Shear  Deformation  Effects 
In  Composite  Laminated  Plates.  (December  1981) 

Brett  Ainsley  Coulter,  B.S.,  Math.,  Colorado  School  of  Mines 

B.S.,  C.E.,  University  of  Colorado 
Chairman  of  Advisory  Comnittee:  Dr.  Y.  Weitsman 

This  thesis  presents  a  closed-form  displacement  solution  for  a 
composite  rectangular  plate  loaded  normal  to  its  surface.  A  high- 
order  plate  theory  with  eight  independent  displacement  functions  is 
used  and  the  mid-plane  deflection  solution  is  compared  with  the 
solution  obtained  from  classical  plate  theory,  which  admits  only 
three  independent  functions. 

The  analysis  provides  the  elastic  solution  for  a  laminated  plate 
made  of  a  graphite/epoxy  material.  Several  laminate  lay-ups,  plate 
aspect  ratios  and  width-to-thickness  ratios  are  considered.  Visco¬ 
elastic  damping,  excluding  inertia  effects,  is  also  considered  with 
frequency  covering  a  wide  spectrum. 
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ABSTRACT 

MODE  I  Delamination  of  a 
Unidirectional  Graphite/Epoxy  Composite 
Under  Complex  Load  Histories  (December  1981) 

James  Scott  Cullen,  B.S.  Civil  Engineering 
University  of  Massachusetts 
Chairman  of  Advisory  Committee:  Dr.  K.L.  Jerina 

A  split  laminate  beam  specimen  is  used  to  obtain  the  relationship 
between  energy  release  rate  and  crack  speed  in  a  unidirectional  graphite/ 
epoxy  composite  under  ambient  and  elevated  temperature  and  humidity  en¬ 
vironments  for  pure  Mode  I  delamination.  The  analysis  employs  a  linear 
elastic  fracture  mechanics  approach  coupled  with  non-linear  beam  theory 
to  account  for  the  large  deflections  and  rotations  produced  experiment¬ 
ally. 
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ABSTRACT 

Compression  Induced  Delamination  in  a 
Unidirectional  Graphite/Epoxy  Composite  (Dec.  1981) 

John  William  Earley,  B.S.  Aeronautical  Engineering 
California  Polytechnic  State  University 
Chairman  of  Advisory  Committee:  Dr.  K.L.  Jerina 

The  opening  mode  delamination  fracture  toughness  of  a  graphite/ 
epoxy  composite  subjected  to  axial  compression  and  a  constant  side 
load  was  investigated  experimentally.  Energy  release  rate  for  a 
stably  growing  crack  (G )  was  determined  using  linear  beam-column 
theory  and  linear  elastic  fracture  mechanics  analysis. 

The  experimental  program  consisted  of  a  symmetric  graphite/ 
epoxy  beam-column  arrangement  simultaneously  subjected  to  axial  and 
lateral  loads.  The  results  from  the  study  indicated  good  agreement 
between  theoretical  and  experimental  displacements,  and  a  good  pre¬ 
diction  of  Mode  I  critical  strain  energy  release  rate  (GjC).  An 
energy  release  rate  of  137.30  N-m/m^  {..784  in-lb/in^)  was  determined 
in  general  for  the  AS! -3502  graphite/epoxy. 
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ABSTRACT 

Mode  I  -  Mode  II  Delamination  Fracture  Toughness  of  a 
Unidirectional  Graphite/ Epoxy  Composite.  (December  1981) 

Peter  Stephen  VanderKley,  B.S.,  University  of  Michigan 
Chairman  of  Advisory  Committee:  Dr.  Walter  L.  Bradley 

The  delamination  fracture  toughness  of  a  graphite/epoxy  composite 
under  varied  mode  I  (opening)/mode  II  (in-plane  shear)  load  ratios  is 
investigated  experimentally.  The  critical  energy  release  rate  G^  jj^c 

is  determined  using  a  double  cantilevered  beam  specimen,  linear  elastic 
fracture  mechanics,  linear  beam  theory,  and  analysis  of  an  asymmetric 
loading  scheme  by  superposition  of  solutions  for  pure  opening  and  pure 
bending  components.  G^j  jj^c  was  found  to  range  from  190  to  590  N/m 
with  the  higher  values  observed  for  greater  Mode  II  contribution.  A 
history  effect  was  also  observed  with  high  mode  II  components  early  in 
the  test  giving  higher  G^j  j^c  values  later  in  the  test.  Fractographs 
of  fracture  surface  indicate  a  correlation  between  fiber/fiber  bundle 
pullout  and  G^j  n)c-  Delamination  cracking  under  mode  II  loads  tra¬ 
versed  entire  ply  thickness  resulting  in  considerable  fiber  breakage 
and  pullout,  which  accounts  for  the  dramatic  increase  in  initial  energy 
release  rate  as  compared  to  pure  mode  I  delamination  where  Gjc  is  found 
to  be  140-160  N/m. 
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ABSTRACT 

Mode  I  Transverse  Cracking  in  an  Epoxy  and  a 
Graphite  Fiber  Reinforced  Epoxy.  (December  1981) 

David  Robert  Williams,  B.S.,  Rice  University 
Chairman  of  Advisory  Committee:  Dr.  Walter  L.  Bradley 

Opening  mode  transverse  fracture  is  studied  in  an 
epoxy  resign  with  and  without  graphite  fiber  reinforcement. 
Energy  absorption  rate  and  critical  stress  intensity  for 
fracture  are  determined  using  compact  tension  specimens 
tested  under  intermittent  fixed  grip  conditions.  The 
graphite/epoxy  composite  manifests  a  significant  increase 
in  fracture  toughness  with  crack  length  as  the  crack  length 
to  specimen  width  ratio  (a/w)  exceeds  0.5.  This  increase 
in  toughness  is  correlated  with  an  increase  in  fiber  break¬ 
age  and  pullout  and  permanent  matrix  deformation.  Results 
obtained  in  this  study  are  compared  with  results  from 
studies  of  delamination  fracture  toughness  of  the  same 
composite  material  to  show  that  transverse  fracture  tough¬ 
ness  is  greater  than  delamination  fracture  toughness  in  the 
graphite/epoxy  composite,  which  is  in  turn  greater  than 
the  neat  resin  toughness  for  the  material  studied. 


8.  PROFESSIONAL  PERSONNEL  INFORMATION 


8.1  Faculty  Research  Assignments 

Each  participating  faculty  member  is  responsible  for  the  research 
conducted  in  at  least  one  specific  area  of  investigation,  as  shown  below. 
In  addition,  most  serve  as  chairmen  of  one  or  more  of  the  graduate  advis¬ 
ory  committees  for  M.S.  and  Ph.D.  students  and,  as  such,  direct  their 
students'  research  project.  The  faculty  also  contribute  to  other  research 
activities  on  the  project  by  serving  on  student  advisory  committees, 
through  technical  meetings,  informal  discussions,  and,  in  some  cases, 
through  specific  research  work. 

The  Principal  Investigator  (R.  A.  Schapery)  has  responsibility  for 
overall  technical  direction  and  coordination  and  for  project  management. 

In  addition  he  has  direct  responsibility  for  certain  research  work,  as 
noted  below. 


Faculty  Member/Departmental  Affiliation 

Dr.  Walter  Bradley/Mechanical 
Engineering 

Dr.  Walter  Hai si er/ Aerospace 
Engineering 

Dr.  Joe  Ham/Physics 

Dr.  Ken  Jerina/Civil  Engineering 
(Now  at  Washington  Univ.,  St.  Louis) 


Primary  Research  Responsibility 

Delamination  and  Fracture 
Properties. 

Development  of  Finite  Element 
Model s 

Resin  Toughening  and  Aging 

Experimental  Data  Base,  Deforma¬ 
tion  and  Failure  Property  Charac 
terization 


Dr.  Richard  Schapery/ Aerospace  and 
Civil  Engineering 


Principal  Investigator  and 
Theoretical  Models  for  Damage 
Growth  and  Fracture 


Dr.  Jack  Weitsman/Civil  Engineering 


Constitutive  Relations, 
Environmental  Effects 
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8.2  Additional  Professional  Staff 
Mechanics  and  Materials  Center 

Mr.  William  Eue  -  Systems  Analyst 

Mr.  Carl  Fredericksen  -  Electronics  Technician 

Mr.  Bob  Harbert  -  Laboratory  Supervisor 


8.3  Spoken  Papers  and  Lectures  at  Conferences  and  Other  Professional 
Activities  of  The  Faculty  Re1ated~to  Composite  Materials 
(1  January  1 981  -  31  December  1 $81 ) : 

W.L.  Bradley 

Invited  Lectures  and  Conference  Presentations: 

"Failure  Analysis:  Methods  and  Case  Histories,"  Series  of  lectures 
given  in  Mexico,  Jan.  1981. 

"Fracture:  Micromechanisms  and  Macroscopic  Measurements  of  Toughness" 
(one  day  presentation  as  part  of  a  short  course  presented  at  TAMU, 
Jan.  1981). 

"Fracture  Toughness  Studies  of  Gray,  Malleable  and  Nodular  Iron," 

AIME  Annual  Meeting,  Chicago,  Feb.  1981. 

"J-Integral  Fracture  Toughness  Studies  in  Gray  Cast  Iron,"  invited 
lecture  for  Metal  Properties  Council,  Feb.  1981. 

"Fracture  of  Commercial  Cast  Irons,"  invited  lecture  to  American 
Foundrymen's  Society  Annual  Meeting,  St.  Louis,  April  1981. 

"Fracture  Behavior  in  Cast  Irons  -  An  Overview,"  invited  lecture 
at  Lawrence  Livermore  Laboratories,  Aug.  1981. 

J.S.  Ham 


Invited  Lectures: 

"Designing  Resins  for  Toughness,"  General  Tire,  Akron,  May  1981 
and  at  Universal  Energy  Systems,  Dayton,  Aug.  1981. 

"Toughness  of  ATX  Resins,"  AFML,  Dayton,  June  1981. 

"Entanglements  in  a  Polysulfone,"  AFML,  Dayton,  Aug.  1981. 

"Physics  and  Polymer  Research,"  McMurray  College,  and  Abilene 
Christian  University,  Abilene,  Oct.  1981;  American  Association  of 
Physics  Teachers,  Fort  Worth,  Nov.  1981. 
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R.A.  Schapery 

Invited  Lectures  and  Conference  Presentations: 

"Nonlinear  Fracture  Analysis  of  Viscoelastic  Composite  Materials 
Based  on  a  Generalized  J  Integral  Theory, 11  Japan-U.S.  Conference 
on  Composite  Materials,  Tokyo,  Jan.  1981. 

"Dynamic  Fracture  of  Viscoelastic  Media,"  ONR  Workshop  on  Solid 
Propellant,  Washington,  D.C.,  Feb.  1981. 

"J  Integral  Theory  for  Nonlinear  Viscoelastic  Media,"  Univ.  of 
Illinois,  Urbana,  April  1981. 

"Fracture  of  Composite  Materials,"  Phillips  Petroleum,  Bartlesville, 
May  1981. 

"Composite  Materials  for  Structural  Design,"  Seventh  Annual  Mechanics 
of  Composites  Review,  Dayton,  Oct.  1981. 

"On  Viscoelastic  Deformation  and  Failure  Behavior  of  Composite 
Materials  with  Distributed  Flaws,"  ASME  Annual  Winter  Meeting, 
Washington,  D.C.,  Nov.  1981. 


Y.  Weitsman 

Invited  Lectures  and  Conference  Presentations: 

"Residual  Thermal  Stresses  in  Cross-ply  Graphite/Epoxy  Laminates," 

22nd  Joint  AIAA/ASME/ASCE/AHS  meeting  in  Atlanta,  April  1981. 

"Environmental  Effects  in  Polymeric  and  Composite  Materials,"  Bell 
Laboratories,  Atlanta,  April  1981. 

"Optimal  Cool-Down  Paths  and  Residual  Thermal  Stresses  in  Cross-Ply 
Composite  Laminates  and  Adhesive  Joints,"  18th  Annual  Meeting  of  the 
Society  of  Engineering  Science,  Sept.  1981. 

"The  Non-Linear  Thermo-Viscoelastic  Characterization  of  Adhesives," 
Univ.  of  Massachusetts,  Amherst,  Oct.  1981. 

"Stresses  in  Adhesive  Joints  due  to  Moisture  Sorption,"  ASME  Annual 
Winter  Meeting,  Washington,  D.C.,  Nov.  1981. 

"The  Thermoviscoelastic  Characterization  of  Polymeric  Materials," 

Bell  Laboratories,  Atlanta,  Nov.  1981. 

"Environmental  Effects  in  Composites,"  the  Technion,  Haifa,  Dec.  1981. 
Technical  Committee  Membership: 

AIAA  Subcommittee  on  Design  Allowables  for  Composite  Materials. 
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In  addition  to  the  above  activities,  the  faculty  attended  several 
conferences  on  composites,  chaired  technical  sessions,  published  papers 
on  other  projects,  and  worked  as  consultants  to  industry. 
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A  Rapidly  Convergent  Scheme  to  Compute 
Moisture  Profiles  in  Composite  Materials 
Under  Fluctuating  Ambient  Conditions 


Y.  Weitsman 


Texas  A&M  University 
Professor,  Mechanical  &  Materials  Center 
Civil  Engineering  Department 
College  Station,  TX  77843 


(Received  February  12,  1981) 


ABSTRACT 

This  paper  presents  a  highly  efficient  numerical  scheme  to  compute  the 
moisture  distribution  in  composite  materials  and  adhesive  joints  under  time 
varying  ambient  relative  humidities  and  temperatures.  The  moisture  diffu¬ 
sion  is  assumed  to  follow  Fick's  laws.  It  is  shown  that  by  appropriate 
switching  among  the  various  forms  of  the  analytic  solutions,  all  involving  in¬ 
finite  series,  it  is  possible  to  attain  extremely  high  accuracy  by  means  of  a 
small  number  of  terms. 

An  example  is  provided  to  illustrate  the  method. 


BASIC  CONSIDERATIONS 


Consider  a  moisture  sorption  process  that  is  described  by  the  classical  dif¬ 
fusion  laws.  In  the  one  dimensional  case  we  have 


dm  _  d2m 
dt  ~  dxl 


(1) 


to  which  we  must  attach  initial  and  boundary  conditions. 

In  (1)  m  =  m(x,t)  is  moisture  content,  x  is  the  spatial  coordinate,  t  is  time 
and  D  is  the  coefficient  of  moisture  diffusion. 

It  has  been  observed  [1,2]  that  the  equilibrium  moisture  content  depends  on 
the  ambient  relative  humdiity,  and  we  shall  also  assume  that  the  boundary 
conditions  are  determined  by  the  same  quantity.  Furthermore,  the  moisture 
diffusivity  was  found  to  be  most  sensitive  to  tempeature  [3,4],  Several  em¬ 
pirical  relationships  were  proposed,  and  we  shall  employ 
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M  =  Cr° 

(2) 

D(T)  =  D„  exp  (A/Tk  ~A/T) 

(3) 

In  (2)  and  (3)  M  is  the  equilibrium  moisture  content,  r  the  ambient  relative 
humidity,  T  the  temperature,  TR  the  reference  temperature  and  A,  C,  a  the 
material  constants. 

In  accordance  with  previous  analyses  [5,6]  we  can  uncouple  the  process  of 
heat  diffusion  from  all  other  time-dependent  material  processes,  e.g.  moisture- 
diffusion  or  stress-relaxation.  This  simplification  is  justified  because  for  all 
practical  temperature  fluctuations  and  geometrical  dimensions  the  time  re¬ 
quired  to  reach  thermal  equilibrium  is  several  orders  of  magnitude  shorter 
than  the  time-scales  for  moisture  diffusion  or  for  relaxation  response.  Conse¬ 
quently,  we  consider  spatially  uniform  temperature  profiles,  namely  T  =  T(t) 
as  prescribed  by  the  fluctuations  in  ambient  temperature,  when  analyzing  tran¬ 
sient  moisture  diffusion. 


SYMMETRIC  EXPOSURE 

Consider  an  infinite  plate  of  thickness  21 .  Let  —  L  <  x  <  L  and  assume  an 
initial  uniform  moisture  distribution  m0.  When  the  plate  is  exposed  to  an 
elevated  ambient  relative  humidity  the  boundary  moisture  is  given  by  n,  name¬ 
ly  m{x  =  ±  L,  t)  -  n.  Due  to  the  symmetry  of  the  present  problem  it  suffices 
to  analyze  only  the  region  0  <  x  <  L . 

For  constant  n  the  moisture  content  m(x,t )  is  given  by  well  known  expres¬ 
sions  [7,8].  Since  we  aim  at  extending  those  expressions  to  the  case  of  fluc¬ 
tuating  ^(t)  and  temperature  T(l)  we  choose  to  represent  them  in  the  following 
form 


m(x,t)  —  mJo(x,t)  =  M  I(x,t)  (4) 

The  functions  /„(*,<)  and  I(x,t)  take  two  alternate  forms 

(  C(x,/)  1  1  -  C(x,f) 

/„(*.<)  =  j  or  ,  /(*,/)  =  <  or  (5) 

'  1  —  E(x,t)  i  E(x,t ) 

In  (5) 


C(x,f)  =  —2  ^  (—  1)"  •  1  cos{P"X/L)  exp(—  p\  t*)  (6) 
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and 


E(x,t)  = 


+  erfc 


(7) 


with  pn  =  (2 n  —  l)n/2 
and  t*  =  Dt/V 


(8) 


The  complementary  error  function  erfc(z)  decays  rapidly  with  z.  Its  asymp¬ 
totic  value  is  given  by  [9]  erfc  znj('/ ’tz)'1  exp(— z2)  consequently,  for  computa¬ 
tional  precision  of  0(10''*)  —  as  obtains  in  “double  precision”  routines  in 
digital  computers  —  we  can  set  erfc  z  =  0  for  z  >  5.877.  In  the  sequel  we  shall 
designate  this  “cut-off”  number  by  A. 

The  rapid  decay  of  erfc  z  implies  that  series  (7)  converges  rapidly  for  short 
times.  On  the  other  hand  it  is  obvious  that  series  (6)  converges  rapidly  for  long 
times.  To  achieve  computational  efficiency  we  should  therefore  switch  among 
the  two  forms  of  equations  (5). 

Straightforward  arithmetic  yields  that  accuracy  of  0(10'“)  is  maintained  by 
the  following  set  of  rules 


for 

(9a) 

)'<"<&) 

|  use  i  terms  in  series  (7) 

for 

(a) 

(?) 

use  five  terms  in  series  (6) 

(9b) 

for  - Q  _  <  t  <  - 2 -  use  I  terms  in  series  (6)  (9c) 

(2 i  +  1)'  (2/  -  IV 

In  (9a)  and  (9c)  i  =  1,  2,  3,  4.  Also  Q  =  14.93  and  A  =  5.877. 

For  t*  >  Q  the  moisture  distribution  is  uniform  to  within  0(10'“). 

It  follows  from  (9a)  that  we  never  need  more  than  the  four  following  terms 
in  series  (7): 
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It  can  be  noted  that  the  form  of  expressions  (9)  remains  valid  for  any 
desired  accuracy  t,  except  that  A  and  Q  depend  on  i.  Obviously  for  a  lesser  ac¬ 
curacy  we  require  even  fewer  terms  in  (6)  and  (7). 

Consider  now  the  case  of  fluctuating  temperatures,  T  =  T(t).  In  view  of  (3) 
the  diffusivity  D  is  now  time  dependent  and  the  non-dimensional  time  /•  in  (8) 
becomes  a  complicated  function  of  real-time  I.  However,  if  we  consider  t*  = 
DKt/L'  at  the  reference  temperature  T  =  T„  then  in  analogy  with  ther¬ 
moviscoelasticity  [10]  we  can  replace  t*  with  the  reduced  dimensionless  time  4* 
whenever  T  =  T(t)  as  follows 

dr  ; 

4*  =  -  /  exp [A/T„  — A/T(s))ds  (10) 

L  0 

The  moisture  distribution  under  fluctuating  temperatures  is  given  by  (4) 
with  t •  replaced  by  4*.  In  view  of  the  single-valuedness  of  4*  =  4*  (0  it  is 
always  possible  to  convert  the  results  back  to  real  time  t . 

Consider  next  the  case  of  fluctuating  ambient  relative  humidity  r  =  r(i).  By 
equation  (2)  this  implies  p  =  p(r).  When  \i  =  n(t)  and  T  =  T(t)  equation  (4) 
yields,  upon  employment  of  the  superposition  integral 


I 

m(x,t)  -  mJo(x,n  =  /  l\x. 4*  (0  -  4*  (t)]  m  '  (r)dr  (11) 

0 

Equation  (11)  must  of  course  be  evaluated  numerically. 

NON  -SYMMETRIC  EXPOSURE 

Consider  now  an  infinite  plate  of  thickness  L  whose  faces  x  =  0  and  x  =  L 
are  exposed  to  different  relative  humidities  which  fluctuate  independently  of 
each  other.  We  still  assume  that  all  temperature  fluctuations  are  spatially 
uniform  within  the  entire  plate. 

The  solution  to  differing,  but  constant  boundary  conditions  m(o,t)  =a° 
and  m(L,t )  =  with  zero  initial  moisture  m(x,o)  =  o  can  be  expressed  as 
follows 
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W>  =  2 


wL(x,c)  = 


erfc 


erfc 


( 

( 


2/i-t-l  —  jr/L 
2/f* 


2/i  +  1  +  jt/Z- 
2/ /* 


) 

) 


(16) 


In  (12 — 16)  t*  =  Df/L’andg,  =  nn. 

Expressions  (14)  are  available  in  the  literature  [8],  while  (IS)  are  obtained 
by  means  of  a  straightforward  Laplace  transform  and  inversion  method. 

Maximal  efficiency  in  evaluating  Ho  and  HL  is  again  obtained  by  switching 
between  their  alternate  forms  given  in  (13)  and  detailed  in  (14—16),  because 
(14)  is  efficient  for  long  times  and  (IS)  is  advantageous  for  short  times.  For  in¬ 
stance,  for  an  accuracy  of  0(10"‘)  we  never  need  more  than  four  terms  in  each 
of  UB,  UL,  Sa  and  SL  as  listed  in  Table  1  below.  For  a  lesser  accuracy  the 
number  of  terms  is  of  course  smaller. 


Table  I.  Number  of  Ttrms  Required  in  Various  Truncated  Series  to  Attain  Accuracy  of 
*/#-'•>  in  Moisture  Profile  a  =  J,*77,  #J  =  lt/n>  log  e) 


Range 

Largest  Number  of  Terms  in  Each  Series 

Total  Number 

■a 

mm 

W0 

u 

s„ 

SL 

Terms 
in  M2) 

0 

0 

0 

0 

0 

0 

? 

1  2  *  T  l 

<T>  <«!»> 

1 

1 

0 

0 

oi 

0 

4 

1  7  *  ?  ’ 
<?*>  “  ‘(f) 

1 

1 

1 

] 

0 

0 

6 

(?»’•  t*S: 

2 

? 

l 

1 

0 

0 

8 

..  . 

-R.-.«t*«  R. 
(<*t)  i 

0 

0 

0 

0 

i 

1 

a.'! 

( i =4 ,3,2,1 ,0) 

For  fluctuating  boundary  conditions  p°(f)  and  */(/),  and  with  varying 
temperature  7X0  we  employ  the  reduced  time  4*  given  in  (10)  and  a  superposi¬ 
tion  integral  analogous  to  (1 1)  to  get 
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m(x,t)  = 


f  j  Ho  (x,  i‘(f)  -  4*(t)] 

<  l  ^ 


(17) 


dRL(r)  dr 
+  H ,  U,4*(0  -  4*(t))  d~~  J  dr 

THE  NUMERICAL  SCHEME 

To  compute  the  moisture  m(x,t )  we  divide  the  time-span  of  interest  t,  into 
n ,  not  necessarily  equal,  sub-intervals.  These  intervals  A;  =  t,  —  t,  _ ,  (/  =  1,2, 
... n )  with  to  =  0  and  /„  =  tf  should  be  selected  in  a  manner  that  both  the  am¬ 
bient  moistures  p°(t)  and  pL(t)  as  well  as  the  temperature  T{t)  are  represented 
to  within  a  satisfactory  approximation  by  the  “staircase”  functions' 

P°(0  =  m,°,  pV)  =  pf.  T(t)  =  T 
for  (18) 

(/  =  1,  2,  ...n) 

Note  that  in  the  symmetric  case  po  =  mo . 

Denote  g,  =  exp(A /TR  —  A  IT)  then  (10)  yields 

4*  =~jj  ^  <'*  0  =  1.2.  .«)  (19) 

*  =  l 


Thereby 

*  =  y  .  1 

(/  =  0,  1,  —  1,  /  =  1,  2,  ...n) 


(20) 


'Obviously,  only  one  ambient  moisture  jj</)  is  involved  in  the  symmetric  case. 
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The  integrals  (11)  and  (17)  are  now  represented  respectively  by  the  sums 


m(x,t,)  =  m0l0(Ar,4,*)  + 


0*.  — M,  _,)  Kx-fy 


(21) 


and 


yt  [o*;-  m;_  ,)  h0  (x, (•)  +  (m/  -  $  _  ,>  hl (x, (,•))  (22) 


Expressions  (21)  and  (22)  remain  valid  for  any  intermediate  time  t,  where 
/,_!</)<  r,,  provided  we  substitute  the  value  of  T  in  place  of  r  in  (19)  and  (20) 
as  well  as  in  (21)  and  (22). 

Computational  efficiency  is  achieved  by  switching  between  the  two  alter¬ 
nate  forms  given  in  (5)  and  (13)  which  is  accomplished  by  testing  the  ranges  of 
4  *  and  4*  according  to  rules  (9a— 9c)  or  in  Table  1,  respectively.  Obviously  (• 
and  4*  must  replace  r*  in  equations  (9)  and  in  Table  1 . 

2A  NUMERICAL  EXAMPLE  FOR  THE  SYMMETRIC  CASE 

To  illustrate  the  method  we  consider  the  case  of  a  sixteen  ply  5208/T300 
graphite/epoxy  laminate  with  L  =  0.04*.  For  this  material  D„  =  1.5019  x 
10"*  inVmin  and  A  =  6340. 

The  composite  laminate  was  considered  to  be  exposed  to  fluctuating  am¬ 
bient  relative  humidity,  which  is  reflected  as  a  fluctuating  boundary  moisture 

and  to  fluctuating  temperatures. 

Specifically,  we  considered  the  situation  in  which  the  ambient  moisture 
level  switched  between  Vi  V#  and  1  Vo  every  5000  minutes,  while  the  temperature 
varied  between  two  fixed  levels  of  350  °K  and  297  °K  also  every  5000  minutes. 
Computations  were  performed  for  moisture  and  temperature  fluctuating  in 
phase  (case  1)  and  out  of  phase  (case  2).  Case  1  is  shown  by  solid  lines  and  case 
2  is  marked  by  dashed  lines  in  Figure  1 . 

The  results,  exhibited  in  Figure  1 ,  show  the  variation  of  moisture  level  with 
time  at  a  station  located  at  x  =  0.035  \ 

Note  that  sharp  slopes  in  m(x,  t)  vs.  t  occur  during  the  high-temperature 
time  intervals.  Consequently  the  in-phase  case  approaches  the  saturation  level 
of  m( x,  t)  =  1 V*  during  the  “wet”  intervals  while  for  the  out-of-phase  case  the 
moisture  level  at  x  =  0.035'  approaches  0.5V#  during  the  “dry”  intervals. 
The  details  are  shown  by  the  heavy  lines  in  Figure  1 . 
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Figure  I.  Moisture  levels  at  x  =  0.035"  vs.  time  in  a  0.08 "  thick  5208/T300  graphite/epoxy 
laminate  that  is  exposed  symmetrically  to  two  cases  of  fluctuating  ambient  relative-humidity  and 
temperature.  Case  1:  R.H.  In-phase  with  temperature.  Case  2:  R.H.  Out-of-phase  with 

temperature. 


ACKNOWLEDGMENTS 

This  work  was  conducted,  in  part,  under  Contract  F33615-79-C-5517  from 
the  Air  Force  Material  Laboratory  (AFML/AFWAL)  and  in  part  under  Con¬ 
tract  F  49620-78-C-003  from  the  Air  Force  Office  of  Scientific  Research 
(AFOSR).  This  support  is  gratefully  acknowledged. 

REFERENCES 

1.  Shen,  C.  H.  and  Springer.  G.  S.,  "Moisture  Adsorption  and  Desorption  of  Composite 
Materials,''  Journal  Comp.  Mat.,  Vol.  10,  Jan.  1976,  pp.  36—54. 

2.  Delasi,  R.  and  Whiteside,  J.  B„  "Effect  of  Moisture  on  Epoxy  Resins  and  Composites."  In 
Advanced  Composite  Materials  —  Environmental  Effects,  STP  658  (ASTM),  J.  R.  Vinson 
—  Editor,  1978,  pp.  2—20. 

3.  Shirrel,  C.  D.,  "Diffusion  of  Water  Vapor  in  Graphite/Epoxy  Laminates."  In  Advanced 
Composite  Materials  —  Environmental  Effects,  STP  658  (ASTM)  J.  R.  Vinson  —  Editor, 
1978,  pp.  21-42. 

4.  McKague,  E.  L.  and  Halkias,  J.  E.,  Private  Communication. 

5.  Weitsman,  Y.,  “Diffusion  With  Time-Varying  Diffusivity,  With  Application  to  Moisture- 
Sorption  in  Composites,”  Journal  Compt.  Mat.,  Vol.  10,  July  1976,  pp.  193—204. 

6.  Douglass,  D.  A.  and  Weitsman,  Y.,  "Stresses  Due  to  Environmental  Conditioning  of  Cross- 
Ply  Graphite/Epoxy  Laminates,"  Proc.  Third  International  Conference  on  Composite 
Materials  (ICCM3),  Paris,  France.  A.  R.  Bunsell  et.  al..  Editors,  Vol.  I,  pp.  529—542, 
Pergamon  Press  (1980). 


357 


Y.  Weitsman 


7.  Luikov,  A.  V.,  "Analytical  Heat  Diffusion  Theory,”  Academic  Press,  1968,  pp.  97—114. 

8.  Crank,  J.,  “The  Mathematical  Theory  of  Diffusion,”  Oxford  University  Press,  2nd  ed., 
1975,  pp.  47,  104—105. 

9.  Abramowiu,  M.  and  Stegun,  L.  A.,  "Handbook  of  Mathematical  Functions,”  National 
Bureau  of  Standards,  June  1964,  p.  298. 

10.  Morland,  L.  W.  and  Lee,  E.  H.,  "Stress  Analysis  for  Linear  Viscoelastic  Materials  With 
Temperature  Variation,”  Transactions  Society  of  Rheology,  Vol.  IV,  1960,  pp.  233—263. 


358 


MM  3724-81-7 
April  1981 


Optimal  Cooling  of  Cross-Ply  Composite 
Laminates  and  Adhesive  Joints 


by 
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Abstract 

This  paper  concerns  the  optimal  cooling  of  symmetric,  balanced, 
cross-ply  composite  laminates  and  of  adhesive  joints  so  as  to  minimize 
the  residual  thermal  stresses  upon  termination  of  the  cool-down  process. 

The  computations  are  based  upon  a  recently  developed  analytical  sch 
and  employ  up-to-date  data  on  graphite/epoxy  laminae.  The  calculations 
consider  the  thermoviscoelastic  response  of  the  polymeric  resins  and 
incorporate  the  temperature  dependence  of  the  coefficients  of  thermal 
expansion. 

It  is  shown  that  the  viscoelastic  behavior  may  contribute  to  a  sign 
ificant  reduction  of  the  residual  stresses. 


♦Professor,  Mechanics  and  Materials  Center 
♦♦Graduate  Student 
Texas  ASM  University 
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Introduction 


This  paper  presents  results  based  upon  a  recently  developed 
T 1  "!*§ 

optimization  scheme1  J  and  data  on  viscoelastic  response  of  graphite/ 
r  31 

epoxy  composites1-  J  whereby  specific  time-temperature  paths  are  selected 
during  the  cool-down  stage  so  as  to  minimize  the  residual  thermal  stresses 
in  cross-ply  laminates  and  in  adhesively  bonded  alumimum/epoxy  joints. 

The  cross-ply  laminates  under  consideration  consist  of  a  symmetric 
and  balanced  lay-up  of  0°  and  90°  plies.  When  such  laminates  are  cooled 
from  the  presumed  stress-free  cure  temperature  down  to  room  temperature 
the  stiffer  fibers  within  each  ply  restrain  the  transverse  shrinkage 
of  the  adjacent,  perpendicularly  directed,  plies  and  these  mutual  geometric 
constraints  introduce  substantial  residual  stresses  in  the  laminate. 

Since  temperature  causes  thermal  expansion  against  geometric  constraints 
it  can  be  viewed  as  a  stress  inducing  agent.  On  the  other  hand,  temperature 
enhances  and  accelerates  the  process  of  stress  relaxation.  The  optimal 
time-temperature  path  achieves  the  best  interplay  between  the  above  mentioned 
competing  effects  so  as  to  minimize  the  residual  thermal  stress  at  the 
termination  of  the  cool-down  phase. 

Elastic  Analysis 

Consider  a  symmetric,  balanced,  cross-ply  composite  laminate,  consist¬ 
ing  of  equal  number  of  0°  and  90°  plies  exposed  to  fluctuating  ambient 
temperature  T(t).  In  view  of  the  sma ; I  hicknesses  of  tie  composite  laminate 
exployed  in  practice  it  is  permissible  tc  ignore  the  process  of  thermal 
diffusion  and  employ  the  approximatic"  mat  T~T(t)  throughout  the  laminate 

*Numbers  in  brackets  indicate  reference  listed  at  the  end  of  this  paper. 

§See  also  ref.  [2]  for  treatment  of  the  same  problem. 


A  similar  circumstance  holds  also  for  adhesive  joints. 


Denoting  by  subscripts  L  and  T  the  longitudinal  (parallel  to  fibers) 
and  transverse  (perpendicular  to  fibers)  directions,  respectively,  we 
have  the  following  stress-strain  relations  in  each  individual  lamina^'^ 


VaLAT  =  r 


TL 


v.j  aT 

eT_aT^I  =  "  +  r 


(1) 


In  (1)  e  and  o  designate  strain  and  stress,  while  a,  E  and  v  denote 
coefficient  of  thermal  expansion,  modulus  and  Poisson's  ratio,  respectively. 
AT  represents  temperature  difference. 

In  view  of  the  particular  lay-up  under  consideration  we  have 

Oy  ”  *"  ~  a ,  Ty  -  -  -  c 


and  straight  forward  manipulations  yield 


a  -  - 


1+v  T  Ty 


(2) 


For  typical  graphi te/epoxy  laminae  we  have  E(~  ?0Ey  (whereby  vLy~  2CVjt^)  . 


With  vlt=*0.28  we  obtain 


a  =  -rEy(ay-a^)AT 


(3) 


Typically  r«0.90 — 0.93. 

Alternately,  consider  two  adherend  plates-  each  of  thickness  h  -  bonded 
together  by  an  adhesive  layer  of  thickness  a.  L.t  the  joint  undergo  a 
spatially  uniform  temperature  excursion  AT, 


-  2  - 


I  J 


Elementary  calculations  then  yield 


a 


e 


1-v 

e 


1 


1-v 


e 


E,  2h 
A 


(4) 


In  (4)  subscripts  "A"  and  "e"  refer  to  the  adherend  and  adhesive, 
respectively.  For  typical  joints  E^/Ee=50  and  h/a  =  100.  Consequently 
(4)  simplifies  to 


a  a;  — —  (a  -a.)M  (5) 

l-v0  6  A 
e 

Viscoelastic  Formulation 

Data  on  graphite/epoxy  laminas^  and  on  "neat"  epoxy  resins^ 

★ 

indicate  that  E^  and  Eg  are  time  and  temperature  dependent.  The  time- 
dependence  can  be  approximated  by  the  form 

E  =  Eo(t+torq  (6) 

while  the  temperature-dependence  can  be  related  by  means  of  a  shift- 
factor  function  a(T),  with  a{Tp)  =  1,  where  customarily  is  room 
temperature  (293°K).  Consequently  (6}  is  extended  to  read 

r  1  "q 

E-  Eo{v[t/a(,)]j  <7> 

Available  data^’^  indicate  that  a(l)  can  be  expressed  by 

a(T)=exp(-  ^  +  B)  (3) 

Employing  the  superposition  integral,  eqns.  (3)  and  (5)  yield^ 

*These  data  are  corroborated  by  many  other  sources  and  publications 
which  are  not  listed  here. 


3 


t  t 


o(t)  =  ctE 


■/ft 


4  t  * 

o  dr 


In  (9)  and  the  sequel,  a  represents  (a^  -  a^)  or  (ag  -  a^)  and  Eq 

absorbs  the  factors  r  or  1/(1 -v  ) .  respectively.  It  may  be  noted  that 

although,  in  a  strict  sense,  r  =  r(t)  the  effect  of  considering  constant  r 

involves  errors  of  no  more  than  1/2%  in  all  practical  cases. 

r  3i 

Data  on  graphi te/epoxy  L  J  indicate  that  aT  and  vary  significantly 
with  temperature.  This  aspect  of  the  material  response  will  be  considered 
later. 

Optimal  Cooling  Paths 

Consider  now  the  case  of  cooling  from  the  initial  cure  temperature 
Tj  down  to  room  temperature  in  a  given  time  interval  t  .  It  has  been 

shown  ^  that  a  minimal  value  of  o(t)  in  (9)  is  obtained  by  following  an 
optimal  path  T(t)  which  includes  an  abrupt  initial  drop  from  Tjto  TQ. 

This  initial  drop  was  given  by  solving  the  transcendental  equation^ 

T  -  T  -  a(T°]  UC 

To  Ti  ■  nr^y 

In  the  open  interval  0  <  t  <  tf  the  optimal  path  T(t)  was  shown 
to  be  smooth  and  continuous,  where  it  is  governed  by  the  nonlinear  integro- 
differential  equation 


dT  E'(t,V  a  ’  ( T(  t ) ) 
at  ■ '  rrtTtjT  irmryawr) 


In  (10)  and  (11)  primes  indicate  derivatives  with  respect  to  the 
argument.  Also 


E(tr  t2)  *  E 


ds 

FOTsT) 


With  given  £(t),  a ( T )  and  T j ,  the  value  of  Tq  is  determined  by 
(10)  and  the  optimal  temperature  path  can  be  constructed  by  iteration 
through  the  employment  of  (11).  In  most  cases  the  final  prescribed 
temperature  7 ^  is  reached  through  a  second  abrupt  step  at  time  t^.. 

A  particular  simplification  occurs  for  a  "power  law"  response,  as 
expressed  in  (6)  and,  in  addi tion, when  the  shift-factor  function  is  given 
by  (8). 

In  this  case  it  was  shown  [1]  that  the  optimal  path  is  given  by 

T  =  C  ln  K  +  Ct  ^ 1 


where,  in  (12) 


C  =  '(q  f  1)A  and  L  =  K  exp ( CTq) 


The  term  K  is  determined  as  the  root  of  the  transcendental  equation 


Ke 


l 

eBto(K  +  Ctf)q 
(q  +  1 )  A 


d 

q  +  1 


(13) 


Under  the  simplifications  which  occur  in  the  present  circumstances 
the  iterative  process  alluded  to  in  |1]  reduces  to  tin-  sraroh  for  tr.e  root 
K  of  eqn.  (13). 

In  view  of  (8)  eqn,  (10)  yields 


T 


0 


-A 


04) 


With  the  optimal  temperature-path  prescribed  by  (12)  the  stresses 


-  5  - 


d 


are  obtained  Dy  direct  substitution  of  (12)  ir.^o  (9).  After  several 
manipulations  we  get 


where 


i 

o(t)  =  aE0  |  j  0(r,t)dT  +  K(TrT0)  Q(0,t)| 


1  _  1 


Q(  • » t ) 


-9  f/K+Ct  Vq  f/K+CtV  q 

.  to  \WJiW)  -u 


+i 


K+C1. 


(15) 


cr(t^)  is  obtained  by  setting  t  =  t^  in  (15).  At  t  =  t^  we 
get 


oU))  ■  °(tp  *  «E0to‘q(Tf  -  V 


(16) 


Temperature-Dependent  Coefficients  o r  Thermal  Expansion 

Recent  data  on  graphite/epoxy  composites  indicate  that  the  quantity 
(ot^-ct^)  increases  monotonical ly  with  temperature.  Consequently 
(9)  must  be  modified  to  read 


alt)  =••  E 


o 


ds 

apT'  V) 


1 

'*■  ( 1  (  d: 
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In  (17)  ct(  T  ( t ) )  =  ot-j.  { T  ( t ) )  -  o .  'T(t)) 
Consider  the  thermal  strain 


eth(t)  *  a(T(t))  AT(t ) 


(18) 


In  view  of  the  one-to-one  correspondence  between  r.th  and  T  we  have 


a ( T ( t ) )  »  6t-th(t) ) 


Therefore  (17)  can  be  re-written  as  follows: 

t  _  t 


0  L  T 


^  U'th(T))  ^ 


Note  that  (20)  is  analogous  to  (9).  Accordingly,  (10)  and 
(11)  yield  the  following  expressions  for  the  optimal  path 


et.h  '  Eth 


I  6(eth) 


dr(t)  E '  '  ( t ,  tf )  8'  (<-  ,  (t) ) 


_ t_h — 

d  t 


" Lf '  p^th' 


In  (22)  E(t1  ,t2)  =  E 


Eqns.  (21)  and  (22)  can  be  expressed  in  terms  of  a(T)  and 
through  the  reintroduction  of  (18)  and  (19).  We  obtain 


Vh 


°<To>»<V _ 

c.(To)a'(To)-a'(To)an0] 


dT 


a 


In  (24)  a  =  a(T( t) )  and  a=a(T(t)). 

The  solution  for  the  optimal  path  T(t)  can  be  generated  by 
means  of  an  iterative  scheme,  as  noted  in  ref.  [1].  However  for  the 
AS/3502  and  T300/5208  graphite/epoxy  systems  the  data  listed  in  Tables  1 
and  2  below  indicate  that  (23)  and  (2 4)  differ  by  at  most  2%  from  expressions 
(10)  and  (11)  which  assume  constant  a.  Consequently,  we  shall  retain 
(12)  as  our  solution  for  the  optimal  path  even  for  the  variable  a  case 
at  hand. 

We  note  however,  that  although  the  optimal  temperature  paths  remain 
unaffected  by  the  variation  of  a,  the  thermal  stresses  do  differ  significantly. 

In  the  particular  circumstance  that 
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The  stress  immediately  following  the  initial  discontinuity  is 


o( o  )  =  Et  "q  a(T  )(T.-T  ) 
'  OO  0  I  0 


(27) 


After  the  final  discontinuity  we  have 


c(t+f)  -  o(t')  +  E0t0_q 


ct(Tf)(Tf-TI)  -  ad^dp-Tj) 


(28) 


Computations 

Calculations  were  performed  for  two  graphite/epoxy  materials  and 
for  an  epoxy/aluminum  joint.  The  data  t3]»[6]  was  reducec|  to  conform 
with  expressions  (3),  (7),  (8),  and  (25)  as  shown  in  Tables  i  and  2 
below: 


Material 

System 

Eo 

(KN/m2) 

r 

q 

(min. ) 

(cm/cm°K) 

A 

(°K) 

B 

AS/3502 

107.055 

.9 

.00775 

— 

1 

aL  =  -  .217xlO-6 

CO 

!  £ 
l  ^ 

i 

46  81  ! 

>, 

ay  =  24. 73x1 0"6 

! 

O 

Cl 

LU 

(at  T  =  372°K) 

<u 

4-> 

I 

£ 

e> 

T300/52o8 

o 

o 

2 

.93 

.00717 

1 

aL=  -  .593xl0-6 
aT=  25.65xl0"6 

6.792 

43.14 

(at  T  37?°K) 

A1  uminum/  Epoxy 

106-54 

1 

.036 

.5 

a  =  5x10  J 
e 

7.91 

37.05 

1 

-  — —  , 

aA  =  2.4x1 0^ 

_ 

Table  1:  Material  Parameters 


Material 

System 

aL 

(cm/cm/°K) 

aT 

(cm/cm/°K) 

“° 

(cm/cm/ °K) 

BUM 

a 

0 

( cin/cm/ °K) 

■  mm 

AS/3502  1  -.484xl0'6 

.72xlO"9 

9.5xl0~6 

.41x1 0~ 7 

T300/5208 

-.418xl0'6 

-.47xl0'9 

11.9xl0-6 

.37xlO~7 

Table  2:  The  Temperature  Dependence  of  the  Thermal 
Expansion  Coefficients.  a(T)  =  aQ  +  a,T. 


In  our  computations  we  employed  Tj  =  450°K,  TR  =  293°K  and  considered 
cooling  times  tx  =  10,  50,  100,  500,  and  1000  minutes. 

The  integrations  (14)  and  (26)  were  performed  numerically.  For  this 
purpose  the  time  interval  t^  was  divided  into  n  equal  sub-intervals 
At  =  t^/n,  resulting  in  intermediate  times  t^  =  i  At  and  intermediate 
temperatures  T.  evaluated  by  means  of  (12).  It  was  found  that  sufficient 
accuracy  was  obtained  with  n=100. 

In  order  to  assess  the  significance  of  the  variation  of  a  with  T, 

a(T.)  +  a(T  ) 

we  computed  (15)  with  a  =  - 2 - ■ 

The  results  are  shown  in  Figs.  1-9.  In  all  those  figures  temperature 
T  is  in  °K,  stress  is  in  MPa  and  time  is  in  minutes. 

Figs.  1  and  2  exhibit  the  optima  cooling  paths  'l'(t),  with  the 
accompanying  residual  stress,  generated  during  the  cool-down  process  of 
typical  aluminum/epoxy  joints  for  cooling  times  t^  =  50  and  100  minutes, 
respectively.  Note  the  jump  discontinuities  at  t.=o  and  t=tf. 

Analogous  results  are  shown  in  Fig.  3  and  4  for  the  0°/90°  T300/5208 
laminate  and  in  Figs.  5  and  6  for  the  0o/90°  AS/ 3502  laminate. 
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In  all  the  above  figures  the  optimal  temperature  paths  were  determined 
by  (12).  We  note  that  those  paths  "undershoot'  the  level  of  room- temperature 
(293°K),  requiring  an  upward  jump  at  t  =  t^. 

The  stresses  in  Figs.  1  and  2  were  computed  according  to  (15)  and 
(16)  employing  constant  a.  On  the  other  hand  the  curves  o(t)  in  Figs. 

3-6  were  obtained  by  employing  { 26 ) - ( 28 )  considering  variable  a  in 
accordance  with  Table  2. 

The  variation  of  the  residual  thermal  stress  o(t^)  with  the  cooling 
time  tjr  is  shown  in  Figs.  7-9.  Also  shown  are  the  variations  of  the  terminal 
temperature  T^  =  T(t~.)  with  t^.  Note  that  a  logarithmic  scale  is  employed 
for  t^. 

For  comparison  purposes  we  also  exhibit  by  dashed-dotted  lines  in 

Figs.  8  and  9  the  final  residual  stresses  o(tf)  which  result  from  considering 

a(T.)  +  oc(fR)  . 

a  constant  ,  average  value  of  aay  =  - 2 -  instead  of  tne  actual 

temperature  dependent  a(T).  Shown  in  dashed  lines  is  the  elastic  value  of 
o(tt),  which  is  computed  with  a  =  a,  . 

T  a  V 

We  observe  that  by  taking  appropriate  account  of  viscoelastic  effects 
it  is  possible  to  attain  reductions  of  about  25 %  in  the  residual  stresses 
relative  to  the  linear  elastic  predictions. 

Remarks  on  Data  Reduction 

The  "power  law"  representation  of  relaxation  or  creep  data  for  polymeric 
resins  or  fiber-reinforced  composites1,"*  in  the  form  F(t.  -  Et  \  can  be 
surmised  from  the  fact  that  the  plots  of  log  E  vs  log  t  are  represented  by 
straight  lines  [3],[6J.  The  slopes  of  those  lines  determine  the  power  q. 
Obviously,  these  plots  cannot  be  extrapolated  to  time  t  -  o,  lest  E(o)~— 
Furthermore,  creep  or  relaxation  data  are  unreliable  at  very  short  times 


11 


r 


because  of  the  dynamic  effects  that  are  technically  unavoidable  at  the 
early  stages  of  experiments. 

Finally,  it  may  also  be  noted  that  for  E(t)  =  Et  ^  the  optimization 

I 

scheme  (11)  fails  because  it  yields  lim  T  (t) — ►  - 

t  -  o 

To  overcome  the  above  mentioned  artifact  select  a  certain  short  time 
t  =  ts>  at  which  the  data  plots  are  smooth  and  beyond  which  E ( t )  =  Et_t*. 

The  data  can  then  be  extrapolated  back  to  t  =  o  employing  Taylor's  expansion 
and  the  result  fitted  into  the  form  E(t)  =  EQ ( t+t^) -C|  of  eqn.  (6). 

It  turns  out  that  by  letting  Eq  =  E(l+q)  and  tQ  =  t&  we  obtain 
representations  in  the  form  of  (6)  that  match  all  data  points  indisting- 
uishably  from  E(t)  =  Et 
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TEMPERATURE 


STRESS  o  (MPa) 


V 


Fig.  3:  The  Optimal  Cooling  Path  T ( t)  and  the  Accompanying  Thermal 
Stress  c(t)  vs.  Time  t  for  T300/5?08  Symmetric, 

Balanced  0°/90°  Laminates  and  Cooling  Time  tf  =  50  min. 


STRESS  ct  (MPa) 


TEMPERATURE  T  (°K) 


g.  5;  The  Optimal  Cooling  Path  T( t )  and  the  Accompanying  Thermal 

Stress  o(t)  vs.  Time  t  for  AS/3502  Symmetric  Balanced  0°/90n 
Laminates  and  Cooling  Time  t^.  -  50  min. 


li 


STRESS  o  (MPa) 


The  Optimal  Cooling  Path  T{t)  and  the  Accompanying  Thermal 
Stress  c(t)  vs.  Time  t  for  AS/3502  Symmetric,  Balanced  0790 
Laminates  and  Cooling  Time  tf  -  100  min. 
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MICROCOPY  RESOLUTION  TEST  CHARI 

NAT'ONA*.  '■'*  >"V.  *•*[■>  •< >• 


TEMPERATURE  T(tp  (°K> 


The  Variation  of  T(t*)  and  the  Residual  Thermal  Stress  o(tf) 
vs.  log  tf  for  the  T300/5208  0°/90°  Laminates.  (The  Dashed- 
Dotted  Line  Corresponds  to  c(tf)  evaluated  with  a  =  a  . ) 
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